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ABSTRACT 


Effects  of  a  transverse  electric  field  on  a  parallel  flow  diffusion 
flame  in  a  flat  combustion  chamber  were  investigated.  Various  mixtures 
of  propane ,  nitrogen,  and  air  were  introduced  separately  at  the  base  of 
an  experimental  combustion  chamber  and  burned  in  a  diffusion  flame  sheet 
located  between  flat  walls  of  the  chamber  which  served  as  anode  and 
cathode.  The  electrode  walls  were  instrumented  to  measure  the  local 
heat  transfer  rate,  local  current  density,  and  pressure. 

It  was  found  that  the  application  of  a  voltage  difference  across 
the  electrodes  moved  ions  out  of  the  burning  zone  and  resulted  in  a 
current  at  the  electrode  walls.  The  heat  transfer  rate  near  the  base 
of  the  flame  was  considerably  increased  on  the  cathode  and  decreased 
on  the  anode;  at  positions  further  from  the  base  of  the  flame  the  elec¬ 
tric  effect  was  lessened. 

Flame  distortion  was  thought  to  be  caused  by  eleccrically  induced 
gas  motion  derived  from  a  gradient  in  current  density  in  the  flow  direc¬ 
tion  and  by  the  onset  of  an  electrically  induced  flame  flickering.  An 
analysis  of  the  products  of  combustion  revealed  that  the  applied  electric 
field  acted  to  increase  the  amount  of  unburnt  solid  carbon  and  decrease 
the  quantities  of  unburnt  fuel  and  carbon  monoxide. 
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The  area  of  study  known  as  electrofluidmechanics  is  concerned  with 
the  behavior  of  a  fluid  in  an  electric  field  when  the  fluid  is  respon¬ 
sive  in  some  way  to  the  electric  field.  One  such  fluid-field  interaction 
has  as  its  basis  the  force  exerted  on  charged  particles  distributed  with¬ 
in  a  region  of  the  fluid,  which,  because  of  molecular  collisions,  result 
in  forces  on  the  fluid.  In  this  type  of  interaction  the  electric  field 
couples  with  the  fluid  mechanic  system  through  the  inclusion  of  an  elec¬ 
tric  body  force  in  the  equations  of  fluid  motion.  Since  it  transpires 
that  combustion  flames  by  their  nature  are  sources  of  partially  ionized 
gas,  a  mechanism  for  the  flame-field  interaction  is  established  and  pro¬ 
vides  the  basis  for  the  study  in  this  investigation. 

The  physical  situation  considered  is  that  of  a  flat  parallel  flow 
diffusion  flame  of  propane/nitrogen  and  air  burning  at  atmospheric  pres¬ 
sure  in  a  flat  combustion  chamber.  Opposing  walls  of  the  combustion 
chamber  serve  as  anode  and  cathode  by  which  an  electric  field  is  impressed 
transverse  to  the  diffusion  flame  sheet.  The  walls  of  the  chamber  are 
instrumented  to  measure  the  internal  local  heat  transfer  rates  with  elec¬ 
trode  voltage,  gas  mixture,  and  gas  velocity  as  operating  parameters. 

The  electrical  measurements  of  local  current  density,  ion  density,  and 
electric  pressure  are  also  made  to  elucidate  the  ways  by  which  the  field 
interacts  with  the  flame.  Exhaust  temperatures  and  compositions  are 
measured  for  all  values  of  the  operating  parameters  so  that  an  overall 
heat  balance  can  be  set  up  for  the  combustion  chamber.  The  objects  of  the 
investigation  are  first  of  all  to  determine  the  physical  effects  of  the 
electric  field  on  the  flame,  with  a  particular  emphasis  on  the  local  heat 
transfer  rates,  and  secondly,  by  an  analysis  of  experimental  data  and 
theoretical  principles  to  determine  the  mechanisms  responsible  for  the 
physical  effects  observed. 

The  range  of  voltages  encountered  in  this  study  are  low  enough  such 
that  the  processes  of  electrical  arcing,  surface  ionization,  and  second¬ 
ary  ionization  do  not  occur,  and  that  the  field  induced  ion  velocity  is 
small  compared  to  the  random  thermal  velocity  of  molecular  motion.  The 
gas  velocities  are  sufficiently  low  that  compressibility  and  magnetic 
effects  are  not  considered. 
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CHAPTER  2 


SURVEY  OF  PREVIOUS  WORK  IN  THE  AREA  OF 
FLAME-FIELD  INTERACTIONS 

Many  experiments  have  shown  that  certain  flames  respond  to  an 
applied  electric  field;  a  review  of  the  earlier  literature  on  the  sub¬ 
ject  is  given  by  Wilson  (l)t  The  typical  experiment  may  be  represented 
by  a  candle  or  Bunsen  flame  located  in  a  horizontal  electric  field  show¬ 
ing  distortion  of  the  flame  toward  the  negative  electrode.  General 
agreement  has  developed  that  the  major  flame-field  effects  are  caused  by 
a  great  disparity  in  the  size  and  mobility  between  the  positive  and  nega¬ 
tive  charge  carriers  which,  in  an  electric  field,  results  in  net  forces 
on  the  flame  gases  in  the  direction  of  the  drift  of  the  ion  with  the  low¬ 
est  mobility  (2). 

Since  the  antecedent  for  the  flame  response  to  an  electric  field  is 
the  production  of  ion  pairs  within  the  flame,  it  is  net  surprising  that 
much  research  has  been  performed  to  establish,  for  a  given  flame  system, 
the  source  of  the  ionization,  the  identification  of  the  ions,  and  their 
number  density. 

With  respect  to  the  source  of  the  ionization,  two  processes  are  con¬ 
sidered:  nonequilibrium  or  chemi-ionization  resulting  from  the  chemical 
reaction  and  equilibrium  thermal  ionization  as  predicted  from  Saha's 
equation  (3)  which  expresses  the  degree  of  ionization  as  a  function  of 
the  temperature  and  ionization  potential  of  the  molecule.  Many  investi¬ 
gate  '4, 5,6,7)  have  measured  ion  concentrations  in  flames  much  higher 

than  predicted  by  Saha's  equation  and  therefore  concluded  that  chemi- 
ionization  was  important.  However,  in  flames  where  material  with  a  rela¬ 
tively  low  ionization  potential  is  present,  such  as  carbon  particles  or 
metallic  salts,  then  thermal  ionization  may  be  large  enough  to  dominate 
chemi-ionization  (8,9,10). 

Identification  of  flame  ions  and  determination  of  their  number 
density  is  still  a  relatively  current  area  of  research.  Mass  spectrom¬ 
eter  and  electrostatic  probe  studies  (6,11,12)  have  shown  for  certain 
hydrocarbon  flames  that  the  most  numerous  positive  ion  is  H30+,  with 
significant  numbers  of  C3H3+.  Number  densities  of  positive  ions  have 
been  measured  from  10s  to  1014  ions  per  cubic  centimeter  in  the  reaction 
zone  (3,5,13,14,15,16,17,18,19,20).  Negative  ions  0"  and  OH"  have  been 
found  to  be  present  in  quantities  two  orders  of  magnitude  lower  than  the 
positive  ion  (ll)»  which  indicates  that  most  of  the  negative  charges  are 
present  as  electrons,  in  agreement  with  other  investigations  (4,6,7,16); 


*  Numbers  in  parentheses  refer  to  references  in  the  reference  list 
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in  some  cases  electron  attachment  to  molecules  may  also  act  to  produce 
negative  ions  (18).  Considerable  increases  in  the  level  of  ionization 
have  been  achieved  by  the  injection  of  metallic  salts  of  metals  with 
low  ionization  potentials,  such  as  cesium  or  potassium  (1,21)  with  a 
consequent  greater  electric  field  interaction  (22). 

In  the  particular  case  of  hydrocarbon  flames,  carbon  formation 
plays  an  important  role  in  contributing  to  flame  ionization.  Nakamura 
(10)  showed  that  there  is  a  close  relationship  between  ion  concentra¬ 
tion  and  carbon  formation  in  a  diffusion  flame.  Others  (9,23)  have 
similarly  found  that  flames  rich  in  hydrocarbon  fuel  produce  much  more 
ionization  than  lean  flames.  One  mechanism  put  forth  is  that  while  the 
carbon  molecule  has  a  relatively  high  ionization  potential,  the  ioniza¬ 
tion  potential  of  an  aggregate  of  carbon  molecules  would  approach  the 
much  lower  value  of  solid  carbon,  and,  therefore,  thermal  ionization 
would  increase.  The  phenomenon  of  carbon  ionization  has  been  suggested 
as  a  means  of  electrically  controlling  flame  carbon  (2,24). 

Although  Sir  J.  J.  Thompson  suggested  as  early  as  1909  (7)  that 
combustion  processes  could  be  controlled  by  the  application  of  a  magnetic 
or  electric  field,  only  recently  has  research  been  initiated  concerned 
with  practical  effects  of  flame-field  interactions,  in  distinction  to 
purely  diagnostic  or  combustion  kinetics  oriented  research.  Payne  and 
Weinberg(2,25)  in  a  preliminary  investigation  discussed  potential  appli¬ 
cations  to  heat  transfer  control  and  carbon  deposition.  A  later  paper 
(23)  discussed  the  limitations  of  flame-field  interactions  imposed  by 
electrical  breakdown  of  the  gases  in  the  electrode  space  and  by  the  rate 
of  ion  generation  within  the  flame.  It  was  noted  that  the  breakdown 
limitation  could  be  optimized  by  suitable  geometric  configurations,  and 
that  the  ion  limitation  could  be  circumvented  by  seeding  with  metallic 
salts. 


Some  interesting  work  has  been  performed  showing  that  electric 
fields  can  under  some  conditions  support  a  flame  that,  if  the  field  were 
removed,  would  be  extinguished.  Both  blow-off  limits  (26)  and  lean  ex¬ 
tinction  limits  (27,28)  were  found  to  be  increased  by  the  application  of 
an  electric  field.  The  maximum  volumetric  heat  release  rate  was  also 
increased  in  an  opposed-jet  diffusion  flame  system  by  the  application  of 
a  field  perpendicular  to  the  flame  surface  (29). 

Ionization  and  a  consequent  electric  field  interaction  was  demon¬ 
strated  to  exist  for  detonation  waves  (30).  The  detonation  velocity 
could  be  either  increased  or  decreased  by  electric  fields,  or  even 
quenched  (31 >32). 

The  energy  supplied  to  the  flame  through  the  application  of  the 
field  in  the  cases  discussed  above  is  small  compared  to  the  chemical 
energy  release  and  is  not  a  direct  influence  in  the  interaction.  There 
has  been  some  interest,  however,  in  augmenting  the  chemical  energy  re¬ 
lease  in  flames  with  electrical  energy  by  means  of  a  superimposed 
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electric  discharge  (33»3^»35)*  The  goal  is  to  generally  increase  flame 
heat  transfer  and  make  more  efficient  use  of  the  chemical  energy,  or  to 
produce  hotter  flames  than  obtainable  by  combustion  alone.  This  concept 
has  the  operational  problem  that  the  electrical  discharge  tends  to  local¬ 
ize  in  the  flame  rather  than  distribute  over  the  region  of  combustion. 


CHAPTER  3 


ELECTRIC  EQUATIONS  OF  A  DIFFUSION  FLAME 
BETWEEN  FLAT  PLATE  ELECTRODES 

3.1  ONE  DIMENSIONAL  ANALYSIS 

3.1.1  Statement  of  the  Problem  -  The  physical  configuration  consid¬ 
ered  is  that  of  a  flat  diffusion  flame  parallel  to  and  located  between 
flat  plate  electrodes  as  shown  in  Fig.  1.  For  an  approximate  one  di¬ 
mensional  analysis  taking  into  consideration  only  the  variation  in  ion 
density  and  field  strength  across  the  combustion  chamber  in  the  x  di¬ 
rection,  the  variation  of  physical  properties  in  the  y  direction  is 
neglected,  and  the  flame  is  considered  to  be  a  flat  sheet  of  uniform 
thickness. 

For  flames  which  are  affected  by  electric  fields,  there  will  be 
ion  pairs  generated  within  the  reaction  zone  because  of  either  thermal 
or  chemi-ionization.  It  would  be  expected  that  the  ion  generation  rate 
would  be  a  maximum  at  some  point  within  the  reaction  zone  and  decrease 
in  both  directions.  For  the  following  analysis  however,  an  ionization 
zone  thickness,  A,  is  taken  as  the  thickness  over  which  a  constant  ion 
generation  rate,  r,  would  produce  the  same  number  of  total  ion  pairs  as 
actually  produced  in  the  flame. 

When  no  electric  field,  is  applied  to  the  flame,  the  number  density 
n+  and  n_  of  the  positive  and  negative  ions,  respectively,  would  be  such 
a  value  that  the  loss  of  ion  pairs  by  recombination  would  equal  the  pro¬ 
duction  of  ion  pairs  by  the  ion  generation  rate  r.  Then  for  no  applied 
field, 


r  -  an+iiw 


(3.1) 


where  a  is  the  recombination  coefficient. 

When  an  electric  field  is  applied  to  the  ionization  zone,  forces 
are  exerted  on  the  ions,  which  tend  to  move  them  out  of  the  ionization 
zone.  For  the  electrode  configuration  in  Fig.  1  the  positive  ions  that 
are  moved  out  of  the  ionization  zone  would  move  across  the  air  zone  to 
the  right  of  the  flame,  and  the  negative  ions  would  tend  to  move  across 
the  fuel  zone  to  the  left  of  the  flame.  Except  for  a  negligible  amount 
of  ions  which  may  be  present  in  the  air  and  fuel  zones  because  of  random 
thermal  ionization,  the  ions  in  the  air  and  fuel  zones  are  either  al  1 
positive  or  all  negative.  Therefore,  no  recombination  can  occur  in 
these  zones,  and  the  ions  would  drift  across  to  the  electrodes  where 
they  are  measured  as  a  current  density,  J0. 
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It  is  assumed  as  an  approximation  that  the  ion  generation  rate,  r , 
and  the  ionization  zone  thickness,  A,  are  properties  of  the  flame  and 
not  affected  by  the  electric  field.  It  is  clear  then  that  an  upper 
limit  on  the  ion  current  density  would  exist  when  the  field  is  strong 
enough  to  remove  every  ion  pair  from  the  ionization  zone  at  a  rate  suffi¬ 
ciently  high  that  recombination  is  negligible.  At  field  strengths  higher 
than  this  value  no  further  increase  in  ion  current  could  be  expected  as 
long  as  there  were  no  secondary  ionization  occurring  in  the  fuel  or  air 
zones  or  on  the  electrode  surfaces.  The  point  at  which  this  condition 
is  reached  is  termed  the  saturation  point  and  the  associated  quantities, 
saturation  current  density,  Js,  saturation  field  strength  Es,  and  satura¬ 
tion  voltage,  Vs,  are  thus  defined. 

At  field  strengths  lower  than  the  saturation  field  strength,  both 
recombination  and  field  removal  are  active  in  removing  ion  pairs  from 
the  ionization  zone.  In  this  case,  the  ion  current  density,  J,  would 
be  expected  to  be  dependent  on  the  applied  field  strength.  The  electric 
analysis  of  the  flame  between  electrodes  can  be,  therefore,  conveniently 
considered  in  two  regimes,  subsaturation  and  saturation,  and  three  zones, 
the  fuel  zone,  the  air  zone,  and  the  ionization  zone,  as  suggested  by 
Weinberg  (23). 


3.1.2  Subsaturation:  The  Ionization  Zone  -  Throughout  this  analysis, 
the  field  strengths  in  the  gas  are  taken  to  be  low  enough  that  the  ions 
may  be  characterized  as  having  a  constant  mobility,  k.  This  means  that 
macro8copically  the  ions  drift  at  a  velocity 


V+  m  k+E 

(3.2) 

V_  m  k_E 

(3.3) 

for  positive  and  negative  ions,  respectively,  where  k  is  the  ion  mobility 
and  E  is  the  field  strength. 

For  the  one -dimensional  treatment,  the  drift  velocity  is  considered 
high  enough  in  relation  to  the  air  and  fuel  velocities  that  the  ions  move 
across  the  air  and  fuel  zones  with  no  appreciable  motion  in  the  y  direc¬ 
tion.  In  other  words,  it  is  taken  that 


Jy+ 

■  Jy- 

.  •  o 

(3.U) 

The  ion  current  density 

is  then 

J+  ■ 

n+v+ 

for 

positive  ions 

(3.5) 

and 

J_  m 

n_v_ 

for 

negative  ions 

(3.6) 

9 


By  combining  Eqs.  (3.2)  and  (3.3)  into  (3.5)  and  (3.6)  there  results 
the  ion  density: 


n+  »  —— 

+  k+E 

for  positive  ions 

(3.7) 

J. 

n_  •  — 

k_E 

for  negative  ions 

(3.8) 

From  published  results  on  flame  ionization  (4,6,7>l6)  there  is  some 
evidence  that  the  negative  ions  are  composed  largely  of  electrons  and 
hence  have  a  mobility  very  large  in  comparison  to  any  conceivable  posi¬ 
tive  ion.  This  means  that  in  the  presence  of  any  electric  field,  the 
right  side  of  Eq.  (3.8)  is  small  compared  to  Eq.  (3.7)  and  hence  the 
negative  ion  density  is  very  small  compared  to  the  positive  ion  density. 
Consequently,  the  negative  ion  density  is  taken  as  zero  in  regions  where 
an  electric  field  exists. 

A  region  of  the  ionization  zone  can  have  either  a  finite  field 
strength  or  a  zero  field  strength.  Considering  a  region  where  the  field, 
is  finite,  continuity  of  positive  ion  flow  gives 


E  ’  r  -  (3.9) 

But  for  regions  where  there  is  a  field  n.  is  negligible  and  Eq.  (3.9) 
reduces  to 


*£  - 

dx 


r 


(3.10) 


For  convenience,  a  coordinate  system  is  placed  at  the  left  side  of  the 
ionization  zone  in  Fig.  1;  i.e.,  let 

s  •  x  -  L  (3.11) 


Then  Eq.  (3.1°)  is  written  as 


dJ 

ds 


r 


(3.12) 


If  the  boundary  condition  on  Eq.  (3.12)  is  taken  as  J  ■  0  at  s  ■  0,  and 
r  is  considered  constant,  the  solution  to  Eq.  (3.12)  is 

J  -  rs  (3.13) 

Equation  (3.13)  above  gives  a  maximum  current  density  at  the  edge  of  the 
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Ionization  zone  equal  to  rA.  This  current  density,  however,  is  the 
maximum  current  density  obtainable  from  the  flame,  contrary  to  the 
assumption  of  subsaturation  conditions.  It  must,  therefore,  be  assumed 
that  for  subsaturation  one  part  of  the  ionization  zone  has  a  zero  field 
strength,  and  another  part  a  finite  field  strength.  Letting  s  ■  sx  be 
the  value  of  s  that  divides  the  ionization  zone  into  the  two  parts,  then 
for  s  <  sx,  E  ■  0  and  J+  ■  0,  and  for  s  >  sx,  Eq.  (3.12)  applies.  The 
details  of  the  following  analysis  are  carried  out  in  Appendix  1  where 
it  is  shown  that 


si  ■ 


(3.14) 


where  J0  is  the  current  density  at  the  edge  of  the  ionization  zone, 
which  is  also  the  current  density  at  the  electrode  suit  faces . 

For  regions  where  s>  sx,  Maxwell's  relation 


dE  m  JJf 
ds  e 

and  the  boundary  condition  E  •  0  at  s  ■  sx  leads  to 


(3.15) 


(3.16) 


(s  >  sx) 

where  k  without  subscript  will  hereafter  be  taken  as  the  positive  ion 
mobility. 

Equation  (3.16)  gives  the  field  strength  at  any  point  within  the 
ionization  zone.  The  field  strength  at  the  edge  of  the  ionization  zone 
facing  the  cathode,  Ef,  is  obtained  by  setting  s  ■  A  in  Eq.  (3.16). 


Ef 


Jo 


(3.17) 


The  potential  drop  across  the  ionization  zone,  V2,  is  obtained  by 
integrating  Eq.  (3.16)  from  s  »  Si  to  s  =  A. 


V2  - 


1 

2 


(3.18) 


n 


/ 


3-1.3  Subsaturation:  The  Air  Zone  -  In  this  region  no  ionization 
or  recombination  occurs  and  J  *  J0.  From  Appendix  1,  using  Maxwell's 
equation 


dE  a 
dx  e 


(3.19) 


and  the  boundary  condition  E  =  Ef  at  s  »  A,  from  Eq.  (3*16)  there  results 


[2Jq  J-A  > 

-  (x-L  +  -  -  A) 

[te  2JS  'J 


(3-20) 


The  potential  drop  acioss  the  air  zone,  V3,  is  obtained  by  integrating 
Eq.  (3*20)  from  x  «  L  +  A  to  x  s  R.  Therefore 

3.1. 4  Subsaturation:  The  Fuel  Zone  -  The  ions  in  the  fuel  zone  are 
composed  largely  of  electrons  and  as  previously  stated  have  a  very  high 
mobility  and  a  negligible  concentration.  From  Maxwell's  equation 


i  .0 

dx  e 


(3.22) 


and  therefore  the  field  strength  is  constant  across  the  zone.  Since 
E  ■  0  at  the  edge  of  the  zone,  E  a  0  everywhere  in  the  fuel  zone.  The 
potential  drop  across  the  zone,  is  also  zero  as  a  result. 


3.1.5  Subsaturation:  Summary  -  The  total  potential  drop  across  all 
zones  must  equal  the  voltage  applied  to  the  electrodes,  V.  Therefore 

The  ionization  zone  is  very  thin  compared  to  the  electrode  spacing. 
(Calcote  and  King  (36)  measured  an  ionization  zone  of  around  one-half 
a  millimeter  thick.)  So  Eq.  (3.23)  can  be  reduced  to 


2  /2£o  \» 

Hke/ 


(R-L)3'2 


(3.24) 


t 

>  1 


where  small  terms  have  been  dropped 
we  find 


Taking  the  derivative  of  Eq.  (3*24), 


dJQ  9keV 
dV  *  4(R-L)3 


(3.25) 


which  will  be  used  later.  Equation  (3*17)  and  (3.24)  can  be  written  in 
nondimensional  form  as 


and 


3.1.6  Saturation;  The  Ionization  Zone  -  In  the  saturation  regime 
all  ion  pairs  generated  in  the  flame  are  moved  out  of  the  flame  by  elec¬ 
tric  forces  and  measured  as  current  density,  Js,  at  the  electrode  walls. 
The  field  strength  in  the  zone  is  obtained  from  Eq.  (3*15)  as  before; 
however,  in  this  case  the  field  at  the  edge  of  the  ionization  zone  is 
unknown.  From  Appendix  1  the  solution  for  the  field  strength  is 


(A2-s2)r  (3.28) 


where  Ef,  the  field  strength  at  the  edge  of  the  zone  facing  the  cathode, 
is  unknown  at  present  and  will  be  determined  later.  The  potential  drop 
across  the  zone,  VZi  is  obtained  from  integration  of  Eq.  (3.28)  from 
s  a  0  to  s  a  A.  From  Appendix  1, 


V2  = 


(3.29) 


3.1.7  Saturation;  The  Fuel  Zone  -  As  before,  the  field  strength 
is  constant  and  equal  to  the  value  at  the  edge  of  the  ionization  zone 
facing  the  anode.  Equation  (3*28)  at  s  =  0  is 


E  = 


(3.30) 
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which  gives  the  field  strength  at  any  point  in  the  fuel  zone.  The  poten¬ 
tial  drop  across  the  zone,  Vx,  obtained  by  integration  of  Eq.  (3.30) 
across  the  fuel  zone  is 


VX  =  Efa  - 


3.1.8  Saturation;  The  Air  Zone  -  The  boundary  condition  on  Eq.  (3*19) 
is  E  »  Ef  at  x  a  L  +  A.  Then  integrating  Eq.  (3.19) » 


E  =|^  (x-L-A)  +  Ef2j  (3.32) 

Integrating  Eq.  (3.32)  across  the  air  zone  to  find  the  potential  drop  V3, 


T 

V 

kc  | 

•  • 


(3.31) 


v3 


kc 

3JI 


2JS(R-L-A) 

kc 


(3.33) 


3.1.9  Saturation:  Summary  -  The  total,  potential  drop  across  all 
zones  must  equal  the  voltage  applied  to  the  electrodes.  Therefore, 


- 

3/2 

-  Ef3 

+  kc 

3JS 

MlH-h-A)  +  E  2 
kc  1 

(3.3*0 


Again  for  A  <  R  and  A  <  R-L,  the  potential  drop  can  be  simplified  to 


(3.35) 


By  nondimensionalizing  Eq.  (3. 35)  in  the  same  way  as  Eq.  (3*27),  there 
results 


V* 


3  Tp.2  M  L  1 

‘  A  T*  e,2a 

J3/2  e,3 

r_5_] 

23/s  L  ]  R-L 

R-lJ  L  2(R-L). 

1  ‘  2s* 

(H-L)J 

(3.36) 


Figure  2  is  a  plot  of  the  nondimensional  field  strength  at  the  flame 
surface  versus  the  nondimensional  voltage  applied  to  the  electrodes,  for 
the  particular  case  of  L/R-L)  »  1  and  A/(R-L)  *  0.04.  The  general  be¬ 
havior  of  the  flame  field  strength  is  that  it  increases  slowly  with  in¬ 
creasing  electrode  voltage  at  low  voltages,  but  at  higher  voltages  the 
field  strength  increases  more  rapidly. 


3.2  WO  DIMENSIONAL  ANALYSIS 

In  the  most  general  case,  the  flame  ionization  rate  and  the  type 
of  ion  may  be  a  function  of  distance  y  along  the  flame.  Suppose  that 
the  quantity  Jjjk  varies  ii  some  way  along  the  flame.  That  is,  let 


^  =  function  of  y  (3*37) 

where  a  subscript  x  or  y  signifies  the  component  of  that  quantity  in  the 
x  or  y  direction,  respectively.  The  one  dimensional  analysis  of  Section 
3.1  will  now  be  assumed  to  be  valid  at  any  y  if  the  corresponsing  Jx  and 
k  from  Eq.  (3*37)  is  used.  Then 


Jx  =*  n+k  Ex  (3.38) 

is  valid  at  any  point.  Writing  Eq.  (3*38)  as 


Jx 

-j—  -  function  of  y  =  n+Ex  (3.39) 

and  substituting  Eq.  (3.19)  yields 


1  =  function  of  y 

2  dx  - Z - 


(3.40) 


Integrating, 

Ex2  .  (function  of  y)x  +  conEt„t  (3.1,1) 

€ 
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which  show3  that  Ex  is  a  function  of  x  and  y.  Substituting  Eq.  (3.41) 
into  Maxwell's  equation  pertaining  to  irrotationality  of  the  field, 


=  (3.42) 

dx 


=  a  nonzero  function. 


Therefore  Ey  is  nonzero  and  Jy  =»  n+kEy  is  also  nonzero,  contrary  to  the 
assumptions  of  the  one -dimensional  treatment.  It  is  concluded  then  that 
if  the  current  density  or  mobility  k  varies  in  the  y  direction,  a  com¬ 
ponent  of  J  and  E  will  result  in  the  y  direction.  A  component  of  E  in 
the  y  direction  will  result  in  electric  forces  and  induced  motion  in 
that  direction,  which  will  later  be  shown  to  affect  the  location  of  the 
diffusion  flame  in  the  chamber. 


£Ex 

ay 

there  results 


CHAPTER  4 


I 

1 

I 

MODEL  OF  A  DIFFUSION  FLAME  BETWEEN  FLAT  PLATE  ELECTRODES 

1 

4.1  DIFFUSION  FLAME  THEORY  WITH  NO 
“?  ELECTRIC  FORCES 

The  case  considered  first  is  the  diffusion  flame  in  Fig.  1  when  no 
..  voltage  is  applied  across  the  electrodes.  The  electrodes  merely  act  as 

confining  walls,  and  the  flame  can  be  analyzed  by  the  method  need  by 
’•  Burke  and  Schumann  (37),  or  by  modifications  of  this  method  (38). 


The  problem  consists  of  determining  the  position  and  shape  of  the 
flame  within  the  combusion  chamber  when  the  flow  rates  of  gases  cure  known. 
Following  the  simplified  treatment  of  Burke  and  Schumann,  the  following 
assumptions  are  made: 

l)  Velocity  is  constant  and  in  the  y  direction  only; 


1! 

] 

3 

] 
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I 
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2)  The  diffusion  of  the  gases  is  one  dimensional  and 
can  be  described  by  Fick's  law  with  a  constant 
diffusion  coefficient  D; 

3)  The  fuel  and  oxidant  have  the  same  density,  so 
gravity  effects  are  negligible. 

Although  assumption  1  is  subject  to  question  considering  that  vis¬ 
cous  effects  tend  to  reduce  the  velocity  near  the  walls  and  that  the 
high  temperature  generated  by  the  flame  tends  to  cause  expansion  and 
higher  velocities  near  the  flame,  it  is  frequently  mentioned  that  assump¬ 
tion  2  tends  to  counterbalance  the  inaccuracies  of  assumption  1  because 
the  diffusion  coefficient  increases  with  temperature  and  tends  to  shorten 
the  flame  while  the  increased  velocity  near  the  flame  tends  to  lengthen 
the  flame.  In  any  case,  the  simple  diffusion  flame  theory  is  discussed 
here  only  as  a  basis  for  the  discussion  of  a  mechanism  for  an  electric 
effect  to  be  discussed  later. 

The  physical  situation  involved  is  the  interdiffusion  of  fuel  and 
air  which  are  initially  separate  at  y  =  0  and  interdiffuse  as  the  mix¬ 
ture  flows  up  the  combustion  chamber  at  constant  velocity  U.  The  fuel 
and  air  are  considered  separately.  Taking  the  air  component  first,  the 
conservation  of  air  in  a  differential  volume  element  leads  to  the  diffu¬ 
sion  equation 


=  D  dfCA 

dt  dx2 


(4.1) 
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Since 


y  =  Ut 


(4.2) 


and  U  is  taken  as  constant,  Eq.  (4.1)  can  be  written 


^Qa  a  2.  ^?CA 

by  U  bxz 


(4.3) 


which  is  the  diffusion  equation  for  the  air  concentration  anywhere  in 
the  combustion  chamber.  The  boundary  conditions  are 


at  X  a  0, 
at  x  a  R, 

at  y  a  0, 
at  y  a  0, 


Ca  =»  0  from  x  a  0  to  L0 
CA  3  *-Ao  from  x  a  L0  to  R 


Equation  (4.3)  with  the  above  boundary  conditions  can  be  solved  by  the 
method  of  separation  of  variables.  The  solution  is 


Ca  -  Cao  (l-Lo/R)  -  f 

nal 


— sin  2S2  cos  aa  exp 
nn  R  R 


(4.4) 


The  concentration  of  fuel,  Cp.,  is  determined  in  the  same  way,  and 


LoCFo 


r2CFo  4  HK^O 

— sin - -  cos 

nal  nn  R 


SB.  exp  \-  (4.5) 

R  L  HR2  J 


The  concentration  of  fuel  and  air  having  been  found  from  Eqs.  (4.4)  and 
(4.5),  the  flame  surface  is  determined  by  the  values  of  x  and  y  such  that 


iCF  *  CA 


(4.6) 


where  i  is  the  number  of  air  molecules  which  combine  with  one  molecule 
of  fuel. 
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4.2  DIFFUSION  FLAME  IN  AN  ELECTRIC  FIELD 
ONE  DIMENSIONAL  MODEL 

In  the  completely  one  dimensional  treatment  of  a  diffusion  flame, 
the  current  density,  flame  location,  and  electric  field  strength  are  not 
dependent  on  y.  Following  the  results  of  Section  3*1>  there  exists  with¬ 
in  the  gas  containing  the  positive  ion  a  body  force  per  unit  volume 

F  =  n+E  =  £  (4.7) 

The  equations  of  motion  then  permit  a  solution 


Vy  =  U  =  constant 

where  vx  and  vy  are  the  respective  x  and  y  component  of  velocity.  The 
body  force  is  effective  in  producing  a  pressure  gradient 

$E-F  (4.8) 

ox 

As  long  as  the  pressure  generated  is  sma33  enough  to  render  the 
effects  of  pressure  diffusion  or  density  changes  small  compared  to  the 
effects  of  concentration  gradients,  which  is  the  case  here,  the  electric 
effect  is  seen  to  produce  only  a  pressure  gradient  and  have  no  first 
order  effect  on  the  shape  or  location  of  the  flame.  There  is  a  question 
of  whether  this  condition  is  stable  and  is  taken  up  in  a  later  section. 


4.3  DIFFUSION  FLAME  IN  AN  ELECTRIC  FIELD 
WO  DIMENSIONAL  MODEL 


The  two  dimensional  treatment  considers  that  the  current  density  J 
obtained  from  the  flame  is  a  function  of  distance  y  along  the  flame.  It 
was  shewn  in  Section  3.2  that  a  variation  in  J  along  the  flame  leads  to 
a  field  strength,  Ey,  and  a  body  force,  Fy,in  that  direction. 

It  is  now  assumed  that  the  impressed  electric  field  strength  in  the 
x  direction  is  large  compared  to  the  induced  component  in  the  y  direction. 
The  body  force  and  the  current  density  in  the  x  direction  are,  therefore, 
large  in  comparison  to  the  components  in  the  y  direction.  From  Eq.  (4.7) 
the  body  force  in  the  x  direction  is 

Fx  .  i  (4.7) 

The  equation  of  motion  is  then 
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P  fw  +  u  |£l=  .  |E  +  £  (4.8) 

L  dx  dy  J  ox  k 

where  inviscid  incompressible  flow  is  considered  and  w  is  the  induced 
velocity  in  the  x  direction.  The  induced  velocity  is  small  such  that 
the  left  side  of  Eq.  (4.8)  is  negligible  compared  to  either  term  on  the 
right.  By  integrating  Eq.  (4.8)  the  pressure  can,  therefore,  be  written 
approximately  as 


P 


Jx 

k 


+  constant 


(4.9) 


The  reference  pressure  is  taken  as  p  s  0  at  the  flame  sheet  where 
x  =  L,  so  Eq.  (4.9)  is 


P  - 


J 

k 


(x-L) 


(4.10) 


The  two  dimensional  aspect  of  the  problem  is  maintained  in  an 
approximate  sense  by  letting  J  be  a  function  of  y  which  could  be  deter 
mined  by  measuring  the  current  density  along  the  electrode.  From  Eq. 
(4.10)  the  pressure  gradient  in  the  y  direction  is 


&  ,  (4.11) 

dy  k  oy 

The  physical  significance  of  the  term  dj/dy  is  now  clear.  A  pressure 
gradient  in  the  y  direction  arising  as  a  consequence  of  the  gradient  in 
current  density  causes  an  acceleration  of  the  gas  on  one  side  of  the 
flame;  an  induced  transverse  velocity  thereby  results. 

For  steady  inviscid  flow  where  the  velocity  in  the  x  direction  is 
small,  Euler's  equation  in  the  y  direction  reduces  to 

£  +  &  =  0  (4.12) 

2  dy  dy 

Substituting  Eq.  (4.11), 


£  duf  +  (x-L)  dJ  =  0  (4.13) 

2  dy  k  dy 
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For  a  case  where  dj/dy  is  constant,  Eq.  (4.13)  integrates  to 


Pu2  +  (x-L)y  dJ  =  ^ 

2  k  dy  *  2 

The  average  value  of  (x-L)  over  the  air  zone  is  (R-L)/2;  then  on  the 
average  Eq.  (4.14)  is 


£uf  +  (£2)1  .  £”!  (4.1?) 

2  2k  dy  2 

It  is  interesting  to  note  that  Eq.  (4.15)  can  be  obtained  in  another 
way  by  first  assuming  that  the  force  field  is  irrotational.  From  irrota- 
tionality 


where 


Sf, 

dy 


1? 

u 

o 

(4.16) 

ox 

n+Ex 

(4.17) 

n+Ey 

(4.18) 

Substituting  Eqs. 


(4.17)  and  (4.18)  into  (4.16)  results  in 


n+ 


a  0 


(4.19) 


The  quantity  in  the  bracket  in  Eq.  (4,19)  is  zero  because  the  field  is 
irrotational,  so  Eq.  (4.19)  reduces  to 


dn+  dm. 

-° 


(4.20) 


As  an  approximation  Ex  and  n+  are  taken  to  be  the  expressions  developed 
in  Chapter  3  for  the  one  dimensional  analysis,  only  now  J  is  taken  to  be 
a  function  of  y. 

From  Eqs.  (3.19)  and  (3.20),  dn+/dy,  dn+/dx  and  Ex  are  calculated 
and  substituted  into  Eq.  (4.20).  For  A  <  <  L,  Eq.  (4.°0)  reduces  to 
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Ey  -  -  [ (4.21) 
L  Jkc  J  oy 

The  body  force  from  Eq.  (4.18)  averaged  over  the  air  zone  is  then 

(4.22) 


(*y)avg  -  -  & 


2k 


8y 


?  the  body  force  is  the  gradient  of  a  potential,  Euler's  equation  in 
the  y  direction  reduces  to 


pug  (R-L)y  dJ  pU* 
2  2k  ay  =  ~ 


which  is  Eq.  (4.15). 


,  A  series  of  approximations  will  now  be  made  to  determine  the  order 
ot  magnitude  of  the  induced  transverse  velocity.  Consider  that  the 
entire  gas  flow  initially  with  uniform  velocity  U  at  y  ■  0  retains  its 
identity  by  assuming  a  hypothetical  interface  between  the  air  and  fuel 
streams.  This  is  not  true  for  the  diffusion  flame  since  the  gases 
immediately  interdiffuse,  and  of  course  the  flame  surface  itself  is  not 
generally  a  streamline.  As  on  order  of  magnitude  approximation  however, 
the  analysis  is  continued. 


For  constant  density  the  interface  must  move  toward  the  cathode  as 
the  air  stream  is  accelerated.  If  x  -  R  -  a  is  the  position  of  this 
interface,  the  interface  velocity  w  is  approximately 


w=-u35 


dy 


(4.23) 


and  for  continuity  of  mass 


au  a  U(R-Lo) 


(4.25) 


Combining  Eqs.  (4.15),  (4.23),  tuid  (4.24),  the  transverse  velocity  of 
the  interface  is 


U2(R-Lo); 


2kp 


f  82J 

.1  y _  + 

dj" 

fus  _  (R-*o)y  bjT*2 

L  ay2 

8y_ 

[  kp  ayj 

(4.25) 
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If  the  velocity  w  represented  by  Eq.  (4.25)  is  interpreted  as  a 
transverse  velocity  away  from  the  flame  in  the  gas  on  one  side  of  the 
flame,  then  this  velocity  opposes  the  diffusion  velocity  of  one  species 
toward  the  flame  and,  therefore,  the  flame  would  be  shifted  toward  the 
cathode. 

In  a  real  flame  the  question  remains  whether  the  effect  of  the 
electric  field  is  a  first  order  effect.  Suppose  for  example  that  a 
diffusion  flame  in  the  combustion  chamber  of  Fig.  1  is  considered  where 
the  fuel  mixture  is  23.6$  propane  and  76.4 $  nitrogen,  and  both  fuel  and 
air  flow  at  a  velocity  of  2.8  inches  per  second.  For  no  electric  field 
present,  the  concentration  profile  and  flame  surface  can  be  determined 
by  the  method  of  Section  4.1.  Once  the  flame  surface  is  determined,  the 
diffusion  velocity  of  any  specie  can  be  determined  for  any  position  along 
the  flame.  If  now  an  electric  field  is  applied  and  the  current  density 
distribution  is  measured,  the  electrically  induced  transverse  velocity 
can  be  determined  by  Eq.  (4.2?).  Suppose  that  bj/by  is  constant  and 
equal  to  a  value  of  0.2  microampere  per  square  inch.  It  will  be  seen 
later  that  this  value  is  representative  of  the  flames  investigated,  and 
is  in  fact  equal  to  the  average  value  of  bj/by  over  the  first  three 
inches  of  one  of  the  flames  later  designated  as  flame  2.  Listed  in 
Table  1  are  the  values  of  the  diffusion  velocity  and  induced  velocity 
for  the  particular  flame  mentioned. 


Table  1 

Diffusion  and  Induced  Velocities  for  an 
Electro-Diffusion  Flame 


Height  above 

Diffusion 

Induced 

Base  of  Flame, 

Velocity, 

Velocity, 

inch 

inch/sec 

inch/sec 

1.48 

0.045 

0.14 

2.9 

0.036 

0.14 

3.0 

0.024 

0.14 

The  values  in  the  table  above  show  that  under  the  assumptions  made 
the  electrically  induced  velocity  is  even  greater  than  the  diffusion 
velocity  of  air  and,  therefore,  a  first  order  effect  would  be  produced 
on  the  diffusion  flame.  The  formulation  of  the  problem  of  determining 
the  flame  location  in  the  electric  field  case  is  achieved  by  expressing 
continuity  of  specie  at  any  point  in  the  combustion  chamber,  whence 


be  _  D  b2C  1  d(Cw) 
dy  U  dx 2  U  dx 


2? 


where  time,  t,  is  eliminated  by  y  =  Ut  as  before,  and  C  is  the  concen¬ 
tration  of  fuel  or  air.  It  is  seen  that  the  last  term  in  Eq.  (4.26)  has 
appeared  in  the  simple  diffusion  equation  as  a  result  of  the  electric 
effect.  The  solution  of  Eq.  (4.26)  requires  a  specification  of  the  in¬ 
duced  velocity  which,  in  general,  is  a  function  of  x  and  y.  If  w  is 
constant,  then  Eq.  (4.26)  has  constant  coefficients  and  is  solved  by  the 
method  of  separation  of  variables. 

The  main  point  of  the  preceeding  discussion  has  been  to  establish 
that  an  electric  field  can  be  expected  to  produce  first  order  effects  on 
a  diffusion  flame  by  means  of  forces  arising  from  a  gradient  in  current 
density  along  the  flame. 
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CHAPTER  5 

THE  QUESTION  OF  INSTABILITY  OF  THE  FLAME 


Chapter  3  considered  a  planar  diffusion  flame  between  parallel  flat 
plate  electrodes  and  the  distribution  of  ion  density,  electric  field 
strength,  and  electric  potential  were  studied.  Since  only  electric  param¬ 
eters  were  investigated,  it  was  assumed  for  simplification  in  that  analy¬ 
sis  that  the  flame  sheet  was  planar  and  stabilized  in  space.  It  is  known 
hower  (39) >  that  under  certain  conditions  electric  forces  can  act  to  de¬ 
stabilize  the  primary  flow  in  fluid  motion,  in  analogy  with  several 
classical  stability  problems  in  fluid  mechanics.  The  question  considered 
at  this  time  is  the  following:  "Given  a  planar  diffusion  flame  Bheet 
stressed  by  a  transverse  electric  field  as  in  Fig.  1,  can  instability  of 
the  flame  sheet  be  expected,  and  what  parameters  serve  as  the  criterion 
for  instability?" 

When  the  flame  sheet  is  in  equilibrium  at  a  certain  distance  from 
the  anode,  there  are  electrical  forces  on  the  flame  sheet  proportional 
to  Ef,  the  field  strength  at  the  edge  of  the  flame.  To  consider  the 
stability,  a  perturbation  in  flame  position  at  constant  electrode  volt¬ 
age  such  as  to  increase  L  is  assumed;  then  if  the  field  strength  Ef  in¬ 
creases  after  the  perturbation,  the  destabilizing  electric  forces  will 
be  increased  and  a  possible  condition  for  instability  exists.  Mathe¬ 
matically,  if  dEf/dL  is  positive,  a  condition  for  electric  destabilization 
of  the  diffusion  flame  sheet  would  exist. 

Since  two  sets  of  electric  equations  exist,  one  for  subsaturation 
and  one  for  saturation,  the  stability  criterion  is  found  for  both  sets. 

Using  Eq.  (3.24)  and  substituting  Eq.  (3.17)  results  in  an  expres¬ 
sion  for  V  in  terms  of  Ef  for  the  subsaturation  case. 


t 


2  r  ^f2js  l 

=  3|_Ake  J 


(R-L)^2 


(5.1) 


( subsaturation) 


The  equivalent  equation  for  saturation  is  Eq.  (3.35).  Taking  d/dL  at 
constant  V  of  Eq.  (5.1)  euid  (3»35)>  we  get  after  resubstituting  Eq.  (3.17) 


I 


dEf 

dL 


(sub saturation) 


(5.2) 
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'A 


and. 


dEf 

P  2Js(R-L)  ‘ 

Ef  +  - 

ke 

'A 

'  2  JSA  ' 
Ef  - 

ke 

‘A 

dL 

EfL 

Ef2- 

“,<4  keEf 

2JS(R-L)  ^ 

Ef2j 

1/4  keEf2 

K€  J 

+  ~ 

L  *£ 

Js 

(saturation) 


(5.3) 


Ef  in  Eq.  (5. 3)  could  be  eliminated  by  using  Eq.  (3.35)  which  gives  Ef 
implicitly  in  terms  of  V  and  Jfl.  Noting  that  Ef  >  JgA/ke  in  saturation, 
it  is  seen  that  dEf/dL  in  Eqs.  (5*2)  and  (5.3)  are  always  positive,  con¬ 
sequently  a  possible  condition  for  instability  exists. 

Consider  now  the  flame  sheet  in  an  equilibrium  position  between 
electrodes  with  a  field  applied  as  in  Fig.  1.  If  the  flame  undergoes  a 
small  displacement,  6,  in  the  positive  x  direction  while  maintaining  its 
equilibrium  ion  density,  a  destabilizing  force,  n+  (dEf/dL)  6,  tends  to 
increase  the  perturbation.  Let  it  be  also  assumed  that  during  the  per¬ 
turbation  the  equilibrium!  profiles  of  air  and  fuel  concentrations  are  not 
disturbed.  Then,  since  the  equilibrium  position  of  the  flame  is  speci¬ 
fied  as  the  surface  of  stoichiometricity,  the  flame  after  perturbation 
will  find  itself  in  a  region  where 


CA  >  iCF  (5.4) 

The  degree  of  departure  from  the  equilibrium  condition  is  proportional 
to  the  quantity 


6  sr 


The  flame  stability  is  governed  by  the  stabilizing  influence  of  the  term 
above  and  the  destabilizing  influence  of  the  term  n+(dEf/dL).  The  stabil¬ 
ity  is  therefore  related  to  the  ratio  of  these  terms,  G,  where 


28 


] 

I 

I 


G 


bEr 

3T 


S  IC*  '  i0F> 


Substituting 


results  in 


n+ 


J 

kEf 


G 


J  dEf 

W^.  HIT 

[CA-iCF] 


(5.5) 


(5.6) 


(5.7) 


The  numerator  of  Eq.  (5*7)  is  different  for  subsaturation  or  satura¬ 
tion.  For  sub saturation,  using  Eqs.  (5.2),  (3.26),  and  (3.27)»  the  num¬ 
erator  becomes 


3JsV*2 
kEf  dL  "  k(R-L) 

(sub saturation) 


and  for  saturation  the  result 


is,  using  Eqs. 


(5.3)  and  (3.26), 


( 


\ 


J  dEf  _  Js 
kEf  dL  k(R-L) 


1  +  -  [l  -  _L> 

A  E'2  E»2 

U  J 

.  —  —  -  ■  x.  ■ ..  ■  , 

-si-  r -i-T  >[— + — T  -  — 

(R-L)  lE*  2-lJ  L(R-L)  (R-L)2J  (R-L) 


> 

/ 


(5.9) 


Where  E'  appears  in  Eq.  (5.9) >  it  may  in  principle  be  eliminated  in  favor 
of  V'  by  Eq.  (3*36)  which  gives  a  relation  between  the  nondimensional 
field  strength,  E',  and  the  nondimensional  voltage,  V',  or  by  use  of 
Fig.  2. 


It  is  interesting  to  investigate  how  the  electric  destabilizing 
parameter  given  by  Eqs.  (5*8)  and  (5.9)  varies  with  applied  voltage. 
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Figure  3  is  a  plot  of  the  nondimensionalized  electric  destabilizing  param¬ 
eters  from  Eqs.  (5.8)  and  (5.9)  for  the  particular  case  where  L/(R-L)  ■  1 
and  A/(R-L)  *  0.04.  In  the  subsaturation  regime  the  parameter  increases 
with  the  square  of  the  voltage.  The  curve  for  saturation,  from  Eq.  (5.9) » 
drops  to  zero  very  near  the  saturation  point  as  indicated  by  the  curving 
down  of  the  curve  near  E*  ■  1  on  Fig.  3*  This  behavior  arises  as  a  result 
of  the  assumptions  leading  to  Eq.  ( 5 •  9)  >  namely,  that  as  a  consequence  of 
taking  the  negative  ion  mobility  as  very  high,  the  potential  drop  across 
the  fuel  zone  is  zero  in  subsaturation,  but  has  a  value,  depending  on  E^, 
above  the  saturation  point. 


The  stabilizing  influence  on  the  flame  sheet  is  the  denominator  of 
Eq.  (5»7)j  related  to  the  concentration  profiles.  In  the  flat  diffusion 
flame  of  the  type  under  consideration,  this  quantity  will  be  a  function 
of  position  along  the  flame.  For  the  flame  used  as  an  example  in  Section 
4.3  (Flame  2  of  the  experimental  flames),  the  stabilizing  parameter  was 
determined  from  a  plot  of  the  concentration  equations  (4.4)  and  (4.5)  and 
is  shown  in  Fig.  4  in  relative  units  as  a  function  of  distance  along  the 
flame.  At  lower  parts  of  the  flame  the  concentration  gradients  are  high 
and  the  stabilization  is  greatest;  at  higher  parts  of  the  flame  the  sta¬ 
bilizing  parameter  decreases. 

Other  models  of  flame  stability  can  also  be  considered;  e.g.,  a 
flame  surface  tension,  oV,  can  be  assumed  where  a  is  the  flame  surface 
charge  density  and  V  is  the  potential  (40).  Stability  of  such  surfaces 
have  been  discussed  in  the  literature  (4l).  Experiments  demonstrating 
the  onset  of  flame  instability  are  discussed  in  a  later  section. 
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CHAPTER  6 

DESCRIPTION  OF  THE  EXPERIMENTAL  COMBUSTION  CHAMBER 


A  parallel  flow  diffusion  flame  combustion  chamber  was  constructed 
to  determine  experimentally  the  mechanisms  and  effects  of  an  electric 
field  applied  transverse  to  the  flow.  The  shape  of  the  cross  section 
in  the  flow  direction  was  rectangular,  with  an  aspect  ratio  of  5»5»  so 
that  a  wide  flat  flame  could  be  produced  within  the  combustion  chamber 
which  would  approximate  the  ideal  flat  diffusion  flame  discussed  earlier. 
The  gaseous  reactants  were  propane  mixed  with  controlled  amounts  of 
nitrogen  as  fuel,  and  air  as  the  oxidant.  The  gases  were  introduced  to 
opposite  sides  of  the  base  of  the  combustion  chamber  as  shown  in  the 
gas  flow  diagram  in  Fig.  5«  As  the  gases  flowed  upward  within  the  com¬ 
bustion  chamber,  the  streams  diffused  together,  producing  a  diffusion 
flame  sheet  when  ignited. 

The  two  walls  facing  the  flame  sheet  were  each  constructed  of  indi¬ 
vidually  cooled  brass  segments,  each  insulated  from  the  adjacent  one  by 
a  l/8-inch  thick  strip  of  asbestos.  Each  brass  segment  was  cooled  by  a 
separate  stream  of  oil  (General  Electric  IOC  transformer  type)  recircu¬ 
lated  by  a  pump.  The  inlet  and  outlet  oil  temperatures  and  flow  rates 
could  be  measured  for  each  segment,  as  well  as  the  temperature  of  the 
segments  themselves.  Both  walls  consisted  of  eight  such  segments,  each 
l/2-inch  thick  and  one  inch  wide.  Above  the  top  segment  of  each  wall 
was  an  additional  l£  inch  wide  segment  which  made  a  total  height  of  10| 
inches  of  cooled  surface  on  each  side.  The  cooled  surface  on  the  fuel 
side  of  the  flame  served  as  the  anode,  and  the  cooled  surface  on  the  air 
side  served  as  the  cathode. 

The  walls  of  the  combusion  chamber  facing  the  edge  of  the  flame  were 
constructed  of  ^-inch  thick  asbestos  board  lined  with  a  l/64-inch  thick 
sheet  of  mica  on  the  edges  adjacent  to  the  electrodes  to  prevent  current 
leakage.  One  edge  wall  had  sections  cut  out  so  that  the  flame  could  be 
observed  in  an  edge  view.  The  inside  of  this  wall  was  fitted  with  a 
sheet  of  l/l6-inch  thick  Pyrex  glass  to  contain  the  combustion  gases. 

It  was  desired  that  the  flame  sheet  be  not  permitted  to  extend  com¬ 
pletely  across  the  inside  width  of  the  chamber  for  two  reasons.  First, 
the  boundary  layer  on  the  glass  window  might  prevent  an  effective  view 
of  the  flame  in  addition  to  the  danger  of  cracking  the  glass,  and  second, 
since  hydrocarbon  gas  was  used,  there  is  the  possibility  of  carbon  deposi¬ 
ting  on  the  edge  walls  and  providing  a  path  for  current  to  flow  between 
the  electrodes.  Consequently,  a  nitrogen  gas  buffer  layer  was  introduced 
at  each  end  of  the  chamber  at  the  base  of  the  fuel  passage.  When  occas¬ 
ionally  carbon  deposits  did  form  on  the  edge  walls,  it  was  noticed  that 
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they  were  more  easily  wiped  off  the  glass  window,  and  for  this  reason 
the  other  edge  wall  was  also  lined  with  glass.  The  cross-sectional 
dimensions  of  the  inside  of  the  chamber  were  l£  inches  between  electrodes 
and  a  width  of  6-7/8  inches,  of  which  5-7/8  inches  could  be  occupied  by 
the  flame  sheet.  A  sketch  of  the  combustion  chamber  configuration  is 
shown  in  Figs.  6  and  7. 

Above  the  electrode  section  of  the  combustion  chamber  was  a  lg-inch 
high  section  of  asbestos  wall,  followed  by  a  one-inch  high  section  into 
which  twelve  iron-constantan  thermocouples  were  inserted  and  spaced  over 
the  cross  section  of  the  chamber  at  that  height.  Above  the  thermocouple 
sides  was  a  sheet  metal  stack  which  tapered  to  a  two-inch  by  two-inch 
cross  section  five  inches  in  height.  In  this  section  a  0.056-inch  diam¬ 
eter  quarcz  tube  was  inserted  to  withdraw  gas  samples  for  composition 
analysif.  Photographs  of  the  chamber  and  instrumentation  are  shown  in 
Figs.  6  and  9. 

In  the  experimental  program  four  different  flames  were  studied;  in 
all  cases  the  various  gas  velocities  at  the  base  of  the  flame  were  uni¬ 
form  for  that  flame.  The  relative  amounts  of  propane  and  air  delivered 
to  the  flame  are  expressed  by  the  equivalence  ratio,  which  is  the  ratio 
of  the  fuel  mass  rate  to  the  fuel  mass  rate  required  for  stoichiometric 
combustion.  The  portion  of  the  air  which  flows  next  to  the  nitrogen 
buffer  layer  is  not  included  in  the  calculation  of  the  equivalence  ratio. 
The  flames  studied  are  listed  in  Table  2. 


Table  2.  Flames  Used  in  the  Experimental  Program 


Equivalence 

Ratio 

Velocity  at 
Base,  in/sec 

Vol  $  Propane 
in  Fuel 

Flame  1 

1.1 

2.8 

18.5 

Flame  2 

1.4 

2.8 

23.6 

Flame  3 

1.7 

2.8 

28.5 

Flame  4 

1.4 

4.2 

23.6 
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Top  View 
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Oil  Passage 


TV  ruo.+i 


A'.rtfc  b  rc.5 


Edge  View 


Fig.  7  Cross  sectional  view  of  the  combustion  chamber 
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CHAPTER  7 


EFFECT  OF  THE  FIELD  ON  THE  APPEARANCE  OF  THE  FLAME 


The  windows  in  the  edge  wall  of  the  combustion  chamber  enabled  an 
observation  of  the  effects  on  the  appearance  of  the  flame  to  be  made. 
Figures  10-13  show  photographs  of  the  flames  taken  through  the  window 
with  a  Polaroid  Model  11 QA  with  close-up  lenses,  ASA  3000  speed  film, 
and  no  special  illumination.  Figures  10,  12,  and  13  were  taken  at  1/125 
second  exposure;  Fig.  11  at  l/60  second  exposure.  Table  3  describes  the 
appearance  of  the  flames  for  no  applied  voltage;  the  trends  among  the 
flames  are  as  expected,  recalling  that  flames  1,  2,  and  3  are  in  order 
of  increasing  fuel  richness,  and  flame  4  is  at  a  higher  velocity. 


Table  3.  Appearance  of  the  Flames  with  No 
Applied  Voltage 


Height  of  Visible 
Flame,  inches 

Height  of  Blue 
Zone,  inches 

Flame  1 

5 

5 

Flame  2 

7 

1* 

Flame  3 

7 

1 

Flame  4 

8 

All  the  photographs  clearly  show  flame  movement  toward  the  cathode 
as  the  electrode  voltage  is  increased.  At  higher  voltages,  flames  2,  3, 
and  4  are  seen  to  impinge  directly  on  the  cathode  a  short  distance  above 
the  base  of  the  flame;  flame  1,  the  leanest  flame,  is  not  as  stabilized 
on  the  dividing  partition  at  the  base  and  is  detached  from  the  partition 
at  an  electrode  voltage  of  2400  volts  (see  Fig.  10,  last  frame).  In 
this  conuition  flame  1  could  be  either  extinguished  by  a  further  increase 
in  voltage  or  reverted  to  its  former  position  by  a  decrease  in  voltage. 

At  low  voltages  where  the  only  effect  was  a  slight  movement  to  the  right, 
no  change  in  luminosity  was  observed;  however,  at  higher  voltages,  when 
the  flame  was  considerably  shortened,  the  apparent  luminosity  was  de¬ 
creased. 

Observation  of  the  flame  as  the  electrode  voltage  was  increased  re¬ 
vealed  a  flickering  of  the  flame  which  was  most  noticeable  above  1600  to 
2000  volts.  The  flame  sheet,  which  previously  had  a  smooth  laminar 
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Fig.  10  Photographs  of  Flame  1  at  Several  Voltages 

0 

0 

0 

y 

u 

42  o 

u 


Fig.  11  Photographs  of  Flame  2  at  Several  Voltages 
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Fig.  12  Photographs  of  Flame  3  at  Several  Voltages 
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appearance  and  was  somewhat  stationary  in  space,  developed  an  unsteady 
appearance.  This  flickering  condition  is  related  to  the  instability  of 
a  flame  in  a  field  and  is  discussed  again  later  in  connection  with  elec¬ 
tric  measurements  on  the  flame. 

It  was  stated  in  Chapter  2  that  carbon  particles  in  flames  may  be¬ 
come  ionized.  Figures  12  and  13  suggest  that  carbon  ionization  is  active 
here;  frames  5  and.  6  show  luminous  parts  of  the  flames  reversing  the 
trend  of  movement  toward  the  cathode  and  instead  moving  toward  the  anode. 
This  leads  one  to  suspect  that  electrons  are  attaching  to  carbon  parti¬ 
cles  and  moving  them  in  the  direction  of  the  anode. 


J 


CHAPTER  8 


ELECTRIC  MEASUREMENTS 


8.1  CURRENT  DENSITY 

Both  the  anode  and  cathode  walls  of  the  combustion  chamber  consisted 
of  brass  segments  which  were  separated  from  each  other  by  a  l/8-inch 
thick  strip  of  asbestos.  Each  segment  was,  therefore,  electrically  in¬ 
sulated  from  adjacent  segments,  and  by  proper  switching  of  the  electrical 
circuit,  shown  in  Fig.  14,  it  was  possible  to  measure  the  current  in  each 
segment  on  both  electrodes  within  0.02  microampere.  Since  only  current 
associated  with  ion  movement  out  of  the  flame  was  desired,  the  possibil¬ 
ity  of  current  leakage  along  the  edge  walls  of  the  chamber  had  to  be 
determined.  Current  measurements  were,  therefore,  recorded  only  when, 
after  turning  off  the  fuel  supply  and  extinguishing  the  flame,  the  current 
dropped  to  a  zero  value.  The  coolant  loop,  which  was  transformer  oil 
circulating  in  plastic  tubing,  was  also  selected  so  that  no  current  leak¬ 
age  would  flow  through  the  coolant  circuit.  Only  current  measurements 
from  the  top  segments  were  not  obtainable  because  of  leakage,  since  the 
upper  part  of  the  chamber  tended  to  collect  carbon  deposits  which  acted 
as  a  current  path,  and  also  the  top  segments  did  not  have  a  conductor 
above  them  at  the  same  potential  to  act  as  a  "guard  ring". 

Four  different  flames  corresponding  to  the  gas  mixtures  and  veloci¬ 
ties  given  in  Table  2  were  formed  in  the  combustion  chamber  and  the 
currents  to  each  segment  of  the  electrodes  were  measured  as  a  function 
of  anode  voltage.  The  cathode,  on  the  air  side  of  the  flame,  was  main¬ 
tained  grounded.  The  current  density  at  the  level  of  a  segment  was  cal¬ 
culated  as  the  current  in  that  segment  divided  by  the  area  of  the  segment 
facing  the  flame.  The  current  densities  thus  obtained  are  presented 
graphically  in  Figs.  15-25. 

One  point  to  be  considered  is  whether,  within  the  range  of  voltage 
applied,  the  saturation  current  was  obtained.  Typically,  it  can  be  seen 
from  the  current  density  figures  that  in  the  upper  parts  of  the  flames, 
the  current  density  first  increased  with  voltage  and  then  leveled  off 
or  decreased,  which  at  first  sight  appears  to  indicate  current  satura¬ 
tion  conditions.  However,  as  shown  in  Chapter  7,  the  shape  and  location 
of  the  flame  is  greatly  changed  at  voltages  above  1600  volts,  and  as  far 
as  upper  parts  of  the  flame  are  concerned,  the  leveling  off  and  decrease 
in  current  density  may  be  due  to  the  fact  that  the  flame  is  shortened 
and  combustion  is  completed  at  a  lower  level  in  the  combustion  chamber. 
Because  of  the  loss  of  the  physical  integrity  of  the  flames  at  high  volt¬ 
age,  it  cannot  be  determined  at  what  voltage  saturation  is  attained. 
Figures  15  and  21,  which  show  the  current  density  near  the  base  of  the 


Fig,  14  Electrical  circuit 
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Electrode  Current  Density,  nicroaape  inch 


Current  Density  Versus  Anode  Voltage,  y  *  5,25  inches 


Fig.  21  Current  Density  Versus  Anode  Voltage,  y  *  0.75  inch 
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Fig.  22  Currant  Density  Versus  Anode  Voltage,  y  =  1*875  Inches 
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Pig.  24  Current  Density  Versus  Anode  Voltage,  y  =  4,125  inches 


flames,  show  no  sign  of  a  saturation  condition,  as  the  current  is  still 
increasing  with  voltage  at  the  highest  voltage  applied.  Considering  the 
flame  as  a  whole,  it  can,  therefore,  be  assumed  that  the  flame  is  sub¬ 
saturated  within  the  limits  of  the  voltages  applied. 

All  of  the  current  density  curves,  Figs.  15-25,  exhibit  the  same 
general  trends.  Near  the  base  of  the  flame,  with  increasing  voltage, 
the  current  increases  slowly  at  first  and  then  shows  a  more  rapid  in¬ 
crease.  At  upper  parts  of  the  flame,  the  current  is  much  lower  and  at 
high  voltages  begins  to  decrease  again  as  previously  mentioned. 

The  current  density  figures  for  flames  1,  2,  and  3»  which  are  all 
at  the  same  velocity,  indicate  that  the  current  density  values  and  their 
changes  with  voltage  show  essentially  the  same  trends.  The  gas  mixture 
then  is  found  not  to  be  an  important  parameter  concerning  flame  current, 
at  least  for  the  range  of  mixtures  and  voltages  used. 

Current  density  Figs.  21-25  show  some  interesting  results  for  flames 
2  and  4,  which  are  the  same  gas  mixture  but  at  different  gas  velocities. 
It  is  seen  that  the  voltage  at  which  the  current  begins  to  increase  or 
decrease  more  rapidly  occurs  at  different  voltages  for  the  two  flames. 

The  trend  is  that  the  higher  velocity  flame  4  "lags 11  in  voltage  behind 
flame  2;  that  is,  the  voltage  where  the  current  shows  a  more  rapid  in¬ 
crease  or  decrease  occurs  at  a  higher  value  for  flame  4.  A  general  con¬ 
clusion  is  that  the  flames  studied  are  affected  more  by  hydrodynamic  and 
electric  variables  than  chemical  variables,  as  may  be  expected  since  we 
are  dealing  with  diffusion  flames. 

The  current-voltage  characteristics  for  the  segments  comprising  the 
cathode  can  be  used  to  obtain  information  concerning  the  mobility  of  the 
positive  ion.  If  it  is  assumed  that  the  negative  charge  carrier  is  an 
electron,  and  the  one -dimensional  equations  of  Chapter  3  are  applicable, 
then  Eq.  (3.25)  applies  to  the  current  density  and  mobility  of  ions  mov¬ 
ing  toward  the  cathode.  Solving  Eq.  (3*25)  for  k,  gives 


k 


(8.1) 


If  the  experimental  current  density  curves  of  Figs.  15-25  are  used  to 
determine  dJ0/dV,  the  mobility  at  various  positions  along  the  flame  may 
be  calculated.  In  this  way  the  mobility  was  found  for  each  flame  at 
various  positions  along  the  flame  using  the  value  of  dJ0/dV  at  800  volts. 
The  results  are  shown  in  Table  4. 
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Table  4.  Positive  Ion  Mobility  Determined  from  Current 
Density  Data,  inch2volt"1sec”1 


Height  Above 
Flame  Base, 
inch 

Flame  1 

Flame  2 

Flame  3 

Flame  4 

0.75 

1.08 

0.20 

0.27 

0.16 

1.87 

0.68 

0.19 

0.53 

0.09 

3.0 

0.31 

0.47 

0.82 

0.28 

4.12 

0.36 

0.42 

O.58 

0.49 

5-25 

0.16 

0.30 

0.21 

0.31 

6.37 

0.09 

0.16 

0.08 

0.11 

7.5 

0.04 

0.14 

0.20 

0.08 

Although  the  mobility  values  in  Table  4  certainly  have  not  been  de¬ 
termined  with  a  great  deal  of  accuracy,  the  results  do  indicate  that  the 
positive  ion  mobility  probably  lies  in  the  range  of  0.1  to  1  inch^olt”1 
sec”1.  As  a  comparison,  the  mobility  values  of  some  singly  charged  posi¬ 
tive  ions  are  (42),  in  the  same  units:  HaO+,  0.17;  C02+>  0.13;  C0+,  0.17; 
N2+,  0.20;  02+,  0.20.  Table  4  does  suggest,  therefore,  that  the  positive 
ion  is  a  singly  ionized  molecule  with  a  molecular  weight  on  the  order  of 
magnitude  of  about  20.  On  the  other  hand,  the  mobility  of  carbon  parti¬ 
cles  found  in  Reference  24  would  be  on  the  order  of  10”3,  and  the  mobility 
qf  the  electron  would  be  of  the  order  of  103  to  104,  in  the  same  units. 

If  the  current  density  data  presented  in  Figs.  12-25  are  plotted  to 
show  the  variation  in  current  density  along  the  electrodes  at  constant 
voltage,  then  an  important  mechanism  for  an  electric  effect  of  the  flame 
comes  to  light.  Figures  26-29  show  the  current  distribution  for  each 
flame  at  1200  and  2000  volts  across  the  electrodes.  For  the  present, 
only  the  distribution  along  the  anode  is  considered.  If  the  major  part 
of  the  negative  ions  are  composed  of  electrons,  as  was  assumed  in 
Chapter  3,  then  because  of  the  large  mobility  a  negligible  coupling 
would  exist  between  the  gas  motion  and  the  electron  current.  It  is  ex¬ 
pected  then  that  the  current  distribution  measured  along  the  anode  re¬ 
flects  the  current  distribution  along  the  flame  sheet. 

It  can  be  seen  from  Figs.  26-29  that  for  all  the  flames  tested  the 
anode  current  is  much  greater  for  the  first  two  inches  or  so  of  flame 
height  than  for  the  rest  of  the  flame  sheet.  As  a  result,  there  is  a 
region  where  the  quantity  dj/dy  has  a  relatively  high  value.  A  large 
gradient  in  current  density  along  the  flame,  dj/dy,  between  the  regions 
of  high  and  low  current  would  give  rise  to  induced  gas  motion  as  dis¬ 
cussed  in  Chapter  4.  In  that  chapter  a  typical  value  of  an  experiment¬ 
ally  determined  gradient  in  current  density  was  used  to  show  that  the 
interaction  with  the  diffusion  process  was  a  first  order  effect. 
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Current  Distribution  along  Electrodes 


Consider  now  the  current  distribution  along  the  cathode,  which 
represents  positive  ion  distribution.  The  most  significant  observation 
from  Figs.  26-29  is  that  the  distribution  of  current  on  the  cathode  is 
much  flatter  than  the  distribution  on  the  anode.  As  previously  mentioned, 
if  the  anode  current  distribution  reflects  the  distribution  at  the  flame, 
then  these  curves  indicate  that  some  positive  drift  downstream  after 
leaving  the  flame  sheet  and  before  reaching  the  cathode. 

Two  effects  may  be  responsible.  First,  the  ion  drift  velocity  may 
be  low  enough  that  the  gas  motion  sweeps  the  ions  downstream;  and  second, 
the  downstream  motion  of  the  ions  may  be  caused  by  an  induced  electric 
field  in  that  direction.  Considering  the  first  effect,  the  average  trans¬ 
verse  electric  field  strength,  Ex,  is  found  by  taking  the  average  value 
of  Eq.  (3.20),  which  results  in,  for  A  «  R-L, 


The  average  drift  velocity  is  kEaVg,  and  the  downstream  movement,  d,  of 
an  average  ion  would  be 


d  - 


where  U  is  the  gas  velocity. 


(8.3) 


Now  the  induced  field  will  be  analyzed  as  to  its  effect  in  displac¬ 
ing  the  positive  ions  downstream.  The  induced  field  given  by  Eq.  (4.22) 
can  be  averaged  over  the  air  zone  to  give 


(By) 


avg  * 


2  f~ 2(R-L)3  T4  c)J 

5  L  Jk  J  by 


(8.4) 


The  downstream  movement,  d,  of  an  average  ion  due  to  this  field  is 


d 


(Ey)avg  (R"^) 
(^x^avg 


3(R-L)g  8j 

5J  by 


(8.5) 


A  comparison  of  the  values  of  the  downstream  movement  of  the  posi¬ 
tive  ions  calculated  from  Eq.  (8.3)  (ion  drag)  and  from  Eq.  (8.5)  (in¬ 
duced  field)  is  given  in  Table  5.  The  experimental  value  was  determined 
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from  Figs.  26-29  by  graphically  finding  the  value  of  y  that  satisfies 


(8.6) 


where  the  second  integral  was  evaluated  over  the  total  length  of  the 
flame.  The  difference  in  the  value  of  y  from  Eq.  (8.6)  for  the  anode 
and  cathode  represents  the  downstream  movement  of  a  ficticious  average 
positive  ion.  The  values  in  Table  5  were  found  assuming  that  the  aver¬ 
age  mobility  was  0.2  inch2volt-1sec"1,  (R-L)  ■  0.75  inch,  and  the  elec¬ 
trode  voltage  was  1200  volts.  The  value  of  dj/dy  was  taken  as  the  change 
in  current  density  between  the  first  and  fourth  segments  of  the  anode 
when  a  voltage  of  1200  volts  was  applied  across  the  chamber. 


Table  5»  Averag*  Downstream  Movement  of  Positive 
Ions,  Inches 


Experimental 

Induced 

Field 

Ion  Drag 

Flame  1 

0.35 

0.14 

0.004 

Flame  2 

2.3 

0.13 

0.004 

Flame  3 

1.6 

0.15 

0.003 

Flame  4 

2.9 

0.35 

0.007 

The  results  in  Table  5  show  that  the  downstream  motion  of  the  posi¬ 
tive  ions  because  of  drag  on  the  flowing  gas  is  much  too  small  to  account 
for  the  experimental  observations.  On  the  other  hand,  the  downstream 
movement  resulting  from  the  calculated  Induced  field  is  much  closer  to 
the  experimental  value.  It  is,  therefore,  concluded  that  these  results 
lend  support  to  the  concept  of  a  field  induced  gas  motion  on  the  air  side 
of  the  flame. 


8.2  ONSET  OF  FLICKERING 

In  Chapter  5  we  discussed  the  conditions  leading  to  instability  of 
the  flame,  and  in  Chapter  7  we  described  observations  of  the  flame  in  an 
electric  field  where  it  was  noted  that  at  certain  voltages  the  flame 
begins  to  flicker. 

A  determination  of  the  onset  of  flickering  was  performed  in  the 
following  way.  A  2.55-million-ohm  resistor  was  placed  in  the  wires 
leading  from  the  first  and  third  segments  of  the  cathode  to  ground. 
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The  current  in  these  segments  is  then  reflected  by  the  voltage  drop 
across  the  resistor;  this  voltage  drop  was  placed  across  the  terminals 
of  an  oscilloscope.  The  voltage  drop  across  the  resistor  was  small 
enough,  on  the  order  of  ten  volts,  that  it  was  assumed  to  not  affect 
the  internal  electric  field  in  any  significant  way. 

For  a  stable,  steady  flame  the  current  would  be  steady  and  no  fluc¬ 
tuation  of  voltage  would  be  sensed  by  the  oscilloscope.  If  flame  insta¬ 
bility  was  generated  by  the  field,  a  fluctuation  in  current  would  be 
expected,  and  a  deflection  observed  on  the  oscilloscope.  Figures  30-33 
show  photographs  of  the  oscilloscope  traces  for  two  positions  of  each 
flame  at  several  values  of  electrode  voltage.  It  is  seen  that  typically 
at  low  voltages  a  straight  line  trace  is  obtained,  indicative  of  stabil¬ 
ity,  and  that  at  higher  voltages  a  definite  fluctuation  is  noted,  indica¬ 
tive  of  flickering.  The  onset  of  flickering  as  indicated  on  these 
figures  is  more  definite  than  observation  with  the  naked  eye.  Figure  3 2 
also  shows  oscilloscope  traces  taken  when  no  flame  existed  in  the  chamber, 
to  prove  by  the  absence  of  fluctuations  that  the  fluctuating  current  was 
due  in  fact  to  the  flame  and  not  associated  with  current  leakage  or  local 
breakdown. 

A  significant  observation  can  be  made  from  the  traces  showing  the 
onset  of  flickering.  The  critical  voltage  at  which  current  fluctuations 
appear  is  lower  for  the  position  at  y  ■  3  inches  than  for  the  position 
at  y  ■  0.75  inches  along  the  flame.  The  approximate  critical  voltages 
are  shown  in  Table  6. 


Table  6.  Voltages  for  the  Onset  of  Current  Fluctuations 


y  ■  0.75  inch 

y  ■  3.0  inch 

Flame  1 

2000 

400 

Flame  2 

2000 

600 

Flame  3 

2000 

800 

Flame  4 

2400 

1200 

It  appears,  therefore,  that  at  some  values  of  voltage,  some  lower  parts 
of  the  flame  may  be  stable  and  ether  parts  at  a  higher  level  destabilized. 
This  agrees  with  Chapter  5  where  it  was  stated  that  steep  concentration 
gradients  at  low  levels  of  the  flame  tend  to  render  those  pants  more 
stable  than  upper  parts  of  the  flame. 

Consider  now  flames  2  and  4,  which  are  the  same  mixture  but  at 
different  velocities.  Table  6  shows  that  the  onset  of  flickering  at 
y  m  3  inches  occurs  at  1200  volts  for  flame  4  and  at  600  volts  for 
flame  2.  A  reason  for  this  difference  is  revealed  by  comparing  the 
values  of  the  stabilizing  parameter  for  the  two  flames.  The  stabilizing 
parameter  is  plotted  in  Fig.  4  for  flame  2  and  has  a  value  of  3»2  at 
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y-  0.?5  inch  y=  3.0  inch 

Flame  1  Flame  1 


Fig.  30  Unset  of  current  fluctuations  in  flame  1 
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y«  0.75  inch 
Flame  2 


y»  3.0  inch 
Flame  2 


Fig.  3i  Unset  of  current  fluctuations  in  flame  2 
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Fig.  32  Onset  of  current  fluctuations  in  flame  3 


71 


y  a  3  inches.  Flame  4  has  the  same  gas  mixture  so  its  curve  correspond¬ 
ing  to  Fig.  4  is  similar;  in  fact,  from  Eqs.  (4.4)  and  (4.5)  we  see  that 
the  concentration  depends  on  the  group  Dy/UR2.  Since  the  velocity  for 
flame  4  is  1.5  times  the  velocity  for  flame  2,  a  height  of  three  inches 
on  flame  2  is  equivalent  to  a  height  of  two  inches  on  flame  4  when  the 
curve  of  Fig.  4  is  used.  The  stabilizing  parameter  of  flame  4  is  then 
found  to  be  about  10$  greater  than  for  flame  2  at  the  same  level.  If 
the  electric  properties  were  the  same,  flame  4  would  consequently  be 
expected  to  require  a  higher  voltage  to  induce  flickering,  which  was  the 
observed  trend. 

A  characteristic  of  the  current-voltage  curves  of  Chapter  8  can  now 
be  explained.  As  the  electrode  voltage  is  increased,  upper  levels  of 
the  *”Lame  begin  to  flicker  first  as  a  result  of  lower  stabilization. 
Accompanying  the  onset  of  flickering  is  a  more  rapid  increase  in  current. 
As  the  voltage  is  further  increased,  lower  parts  of  the  flame  flicker 
and  cause  a  more  rapid  increase  in  current  there.  High  current  density 
gradients  then  cause  induced  gas  motion  which  results  in  combustion 
being  completed  at  a  lower  level  of  the  combustion  chamber;  as  a  result 
the  currents  in  the  upper  parts  of  the  chamber  decrease. 


8.3  ELECTRIC  PRESSURE 

Holes  of  0.05-inch  diameter  were  drilled  through  the  walls  compris¬ 
ing  the  anode  and  cathode  at  a  position  midway  along  the  width  of  the 
combustion  chamber  and  at  a  level  of  2.19  and  5.56  inches  above  the  base 
of  the  flame.  The  holes  in  the  anode  and  cathode  at  the  same  level  were 
directly  opposite  each  other,  and  served  as  static  pressure  taps  to 
measure  the  difference  in  pressure  across  the  flame  as  a  function  of 
electrode  voltage.  The  pressure  measuring  instrument  was  a  Flow  Corpora¬ 
tion  micromanometer  fitted  with  a  special  traversing  microscope  capable 
of  reading  relative  pressures  to  within  1.4  x  10"5  inches  of  water. 
Because  of  zero  shift,  the  pressures  are  only  accurate  to  about  3  x  10”5 
inches  of  water. 

Considering  first  a  region  between  the  flame  and  the  cathode,  in 
the  presence  of  an  electric  field  there  is  a  body  force 

F  -  i  (8.6) 

existing  on  the  gas.  The  body  force  appears  in  the  equation  of  motion 
in  the  x  direction  to  yield 
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where  w  is  the  induced  transverse  velocity,  u  is  the  velocity  parallel 
to  the  electrodes,  and  incompressible,  inviscid  flow  is  considered.  The 
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terms  on  the  left  side  of  Eq.  (8.7)  are  several  orders  of  magnitude  less 
than  the  terms  on  the  right  side  for  the  induced  velocities  considered, 
so  Eq.  (8.7)  is  rewritten  as 


8p,  £ 

8x  k 


(8.8) 


Since  the  pressure  and  current  density  are  measured  quantities,  Eq.  (8.8) 
serves  as  a  method  to  determine  the  mobility.  Figure  34  shows  a  plot  of 
the  cathode  pressure  minus  the  anode  pressure  versus  voltage  for  the 
position  of  2.19  inches  above  the  base  of  the  flame.  These  curves  show 
only  the  electric  pressure  and  therefore  indicate  zero  pressure  differ¬ 
ence  at  zero  electrode  voltage.  Actually,  there  was  a  small  measured 
difference  between  the  electrode  pressures  at  zero  voltage  equal  to  1.55 
for  flame  1,  0.28  for  flame  2,  2.4  for  flame  3>  and  1.0  for  flame  4  where 
the  value  represents  cathode  pressure  in  excess  of  anode  pressure  in 
inches  of  water  times  10"4. 

Since  it  is  assumed  a  priori  that  the  negative  ion  at  this  level  is 
an  electron,  the  pressure  gradient  between  the  anode  and  the  flame  is 
very  small.  The  positive  ion  mobility  from  Eq.  (8.8)  is  then 


k  .  (8.9) 

where  P  is  the  cathode  pressure  in  excess  of  the  anode  pressure  from 
Fig.  3*+*  The  distance  (R-L)  was  not  directly  measured,  but  from  visual 
observations,  ‘v,e  following  are  taken  as  approximate  values:  for  400 
volts,  1  Inch;  .or  800  volts,  7/8  inch;  for  1200  volts,  3/4  inch;  for 
1600  volts,  l/2  inch;  for  2000  volts,  l/4  inch.  Using  the  measured 
values  of  electric  pressure  from  Fig.  34  and  measured  value  of  current 
density,  the  mobility  was  calculated  using  Eq.  (8.9)  and  is  listed  in 
Table  7  below. 


Table  7.  Positive  Ion  Mobility  at  y  ■  2.19  Inches  Determined  From 
Electric  Pressure  Measurements,  inch2volt"1sec"1 


Electrode 

Voltage 

Flame  1 

Flame  2 

Flame  3 

Flame  4 

4oo 

0.80 

0.30 

800 

0.39 

0.27 

0.07 

0.64 

1200 

O.56 

0.22 

O.36 

O.65 

1600 

O.65 

0.32 

0.34 

0.87 

2000 

0.67 

0.41 

0.23 

0.21 

The  values  in  Table  7  are  in  the  same  range  as  the  values  of  mobil¬ 
ity  found  earlier  and  listed  in  Table  4.  There  is  a  strong  indication, 
therefore,  that  the  assumption  of  a  positive  ion  with  an  atomic  weight 
on  the  order  of  about  20  and  a  negative  ion  being  an  electron  is  justi¬ 
fied  for  the  portion  of  the  flame  2.19  inches  above  the  base. 

Attention  is  now  turned  to  the  pressure  measurements  at  a  level  of 
5.56  inches  above  the  base  of  the  flame,  which  are  shown  in  Fig.  35*  It 
is  readily  apparent  that  conditions  have  changed,  for  at  this  level  the 
electric  pressure  on  the  anode  is  greater  than  on  the  cathode— a  reversal 
of  the  effect  observed  at  the  lower  level  of  the  flame.  Apparently  the 
average  negative  ion  at  this  level  is  not  an  electron  but  is  an  ion  with 
a  mobility  even  lower  than  the  positive  ion.  By  writing  the  pressure 
rise  across  the  aiyl  and  fuel  zones  we  get 
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If  the  positive  ion  mobility  is  assumed  to  be  unchanged,  Eq.  (8.10)  leads 
to  values  of  negative  ion  mobility  of  around  0.1  inch2volt"1sec"1,  or  on 
the  order  of  the  value  of  the  positive  ion  mobility.  Of  course,  since 
only  average  mobilities  are  detected,  some  negative  ions  may  have  a 
mobility  larger  than  this,  and  others  may  have  a  mobility  lower. 

Two  mechanisms  for  the  increase  in  negative  ion  mobility  along  the 
flame  may  be  active.  One  case  is  electron  attachment  to  neutral  mole¬ 
cules.  Near  the  base  of  the  flame  the  gas  between  the  anode  and  the 
flame  is  mostly  nitrogen  with  some  propane,  and  small  amounts  of  combus¬ 
tion  products.  Nitrogen  has  a  very  low  probability  of  attachment  and 
does  not  act  as  a  charge  carrier.  The  products  of  combustion,  notably 
water  vapor,  however,  readily  attach  to  electrons  and  can  act  as  an 
electron  carrier.  As  more  combustion  products  are  present  at  higher 
levels  in  the  chamber,  an  increasing  percentage  of  the  negative  ions 
would  have  a  mobility  corresponding  to  ionized  molecules. 

Another  process  for  decreasing  the  mobility  is  the  attachment  of 
electrons  to  carbon  particles.  Carbon  particles  would  be  formed  in  the 
combustion  chamber  by  pyrolysis  of  fuel  in  the  region  near  the  flame. 

The  particles  would  be  in  a  region  of  electron  flux.  Weinberg  and  Place 
(24)  have  shown  in  a  discussion  of  electrical  control  of  flame  carbon 
that  under  the  above  conditions  thermal  ionization  of  carbon  particles 
and  electron  attachment  to  carbon  particles  are  both  active  processes 
which  produce  charged  particles  of  opposite  sign.  Honda  (42)  has  also 
shown  that  in  a  soot-producing  flame,  the  soot  had  an  electric  charge. 
Since  a  net  electric  pressure  was  measured  on  the  anode  at  this  level 
of  the  chamber,  it  is  thought  that  the  electron  attachment  is  predomi¬ 
nant. 
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Anode  Voltage,  volte 


Electric  pressure  at  y  -  5.562  Inches  above  the  flame  base 


This  increase  of  pressure  on  the  anode  agrees  qualitatively  with 
observations  of  the  flames  discussed  in  Chapter  7.  It  was  observed  that 
movement  of  luminous  gases  toward  the  anode  took  place,  reversing  the 
general  trend  of  movement  toward  the  cathode. 

The  possibility  that  pressure  differences  on  the  electrodes  may  be 
due  to  an  unsymmetric  heating  of  the  flowing  gases  has  not  been  estab¬ 
lished,  although  this  effect  may  be  present.  However,  since  the  mobil¬ 
ity  values  calculated  from  the  pressure  differences  at  y  *»  2.19  inches 
gave  expected  values,  it  seems  that  the  pressure  difference  is  caused 
by  current,  at  least  for  that  location. 


8.4  ELFCTROSTATIC  PROBE  MEASUREMENTS 

A  probe  technique  was  used  to  measure  the  local  ion  density  at  a 
point  within  the  combustion  chamber  i:i  the  air  zone  between  the  flame 
and  the  cathode.  A  steel  wire,  0.008-inch  in  diameter,  was  stretched 
horizontally  at  a  level  of  2.19  inches  above  the  base  of  the  flame 
and  parallel  to  the  cathode  wall  at  a  distance  of  3/8  inch  from  the 
cathode.  The  wire  was  supported  at  each  end  by  a  piece  of  quartz  tubing 
projecting  through  holes  drilled  in  the  cathode  wall,  and  cemented  in 
place.  The  wire  supports  were  5-1/8  inches  upert  and  hence  the  wire 
stretched  along  practically  the  entire  width  of  the  flame  sheet. 

The  wire  probe  was  connected  to  a  probe  measuring  circuit,  Fig.  36, 

■  which  the  current  to  the  probe  could  be  measured  as  a  function  of  the 
voltage  impressed  on  the  wire.  From  simplified  probe  theory,  (43),  the 
ion  density  and  space  potential  at  the  position  of  the  probe  can  be  de¬ 
termined  by  the  probe  curren-voltage  characteristic.  For  the  case  of 
positive  ions  and  neutrals  only,  the  current- voltage  characteristic 
plotted  on  a  semi-log  graph  is  typically  as  shown  in  Figs.  37-40.  The 
significant  point  on  such  curves  is  the  point  where  the  theoretically 
linear  portion  of  the  curve  at  high  positive  voltage  breaks  off  to  a 
flatter  curve  at  high  negative  voltages.  At  the  point  corresponding  to 
the  knee  of  the  curves,  the  voltage  is  equal  to  the  space  potential  at 
the  point  of  the  probe,  and  the  probe  neither  attracts  nor  repels  ions 
electrically.  The  probe  current  at  this  condition  is  then  a  measure  of 
the  local  ion  density.  If  Jp  is  the  probe  current  density  at  the  knee 
of  the  curve,  K  is  Boltzmann’s  constant,  and  T  is  the  temperature,  the 
positive  ion  density  is  given  by 

-  -  "MM  f  (8-u) 

The  current-voltage  characteristic  for  the  probe  for  each  flame  at 
a  voltage  of  800  volts  across  the  combustion  chamber  is  shown  in  Figs. 
37-40.  The  probe  currents  fluctuated  during  the  measurement,  and  the 
range  of  fluctuation  was  recorded  and  plotted  as  a  vertical  line.  As 


Fig.  36  Electrostatic  probe  circuit 
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Fig.  40  Electrostatic  Probe  Characteristic  for  Fl«me  4 


usual,  the  knees  of  the  probe  characteristics  are  not  clearly  defined, 
but  a  representative  point  can  be  chosen  such  that  the  value  is  correct 
within  about  an  order  of  magnitude  at  most.  The  ion  density  was  thus 
determined  from  the  curves  in  Figs.  37-^0  and  the  results  shown  in 
Table  8.  An  atomic  weight  of  20  and  a  temperature  of  200°F  was  used  to 
evaluate  Eq.  (8.11). 


Table  8.  Positive  Ion  Density,  Ions  Inch" 


Flame  1 
Flame  2 
Flame  3 
Flame  4 


Probe 

Measurement 

0.09  x  10e 
0.1  x  10s 
0.31  x  109 
0.61  x  10° 


Equation  8.12 

1.7  x  10s 
3.0  x  109 
4.5  x  10® 
5.4  x  108 


The  ion  density  can  also  be  calculated  using  the  one-dimensional 
electric  equations  of  Chapter  3-  Equation  (3.20)  for  A<<L  leads  to 
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Using  a  mobility  value  of  0.2  inch2volt"1sec"1  and  (x-L)  ■  l/2  inch,  and 
measured  values  of  current  density,  the  ion  density  was  calculated  from 
Eq.  (8.12)  and  listed  in  Table  8.  The  probe  data  is  seen  to  indicate  a 
lower  ion  density  than  Eq.  (8.12)  by  about  one  order  of  magnitude.  In 
dealing  with  experimental  measurement  of  number  densities  of  ions  however, 
agreement  within  a  factor  of  10  is  occasionally  accepted  as  satisfactory. 


CHAPTER  9 


EFFECT  OF  THE  FIELD  ON  THE  PRODUCTS  OF  COMBUSTION 


The  question  arises  when  dealing  with  electric  field  effects  on 
flames  whether  the  field  alters  the  combustion  process  itself,  therby 
possibly  changing  the  heat  release,  composition  of  products,  or  tempera¬ 
ture. 


The  mechanisms  by  which  the  applied  field  can  affect  the  combustion 
process  can  be  categorized  into  two  types.  In  one  type  we  can  list  hydro- 
dynamic  and  mechanical  effects  caused  by  electric  forces  on  the  gases  or 
particles,  as, for  example,  the  induced  gas  motion  discussed  in  Chapter  4. 
It  can  be  expected  that  hydrocarbon  diffusion  flames  of  the  type  under 
investigation  would  exhibit  a  particularly  strong  influence  by  mechanical 
effects  because  of  the  existence  of  a  pyrolysis  zone  near  the  flame  sheet 
where  fuel  is  heated  and  chemically  changed.  Induced  gas  motion  would 
conceivably  affect  the  transfer  of  heat  into  the  pyrolysis  zone,  thereby 
directly  changing  the  amount  of  fuel  pyrolyzed,  or  carbon  particles 
could  become  charged  and  consequently  move  through  the  combustion  chamber 
under  the  influence  of  electric  forces.  The  subsequent  growth  or  com¬ 
bustion  of  the  particle  would  then  be  determined  by  electrical  as  well 
as  chemical  and  thermal  conditions.  This  was  demonstrated  in  Reference 
24  where  the  field  strength  was  adjusted  to  hold  a  carbon  particle  stat¬ 
ionary  within  the  pyrolysis  zone,  resulting  in  large  visible  growth  of 
the  particle. 

The  other  type  of  interaction  with  the  combustion  process  is  the 
effect  on  chemical  kinetics  during  combustion.  Ionized  intermediates 
and  radicals  formed  during  the  complex  chain  reactions  could  be  dis¬ 
placed  in  such  a  way  as  to  a-  ter  the  final  reaction.  Unstable  reaction 
species  could  be  moved  to  lower  temperature  regions  where  nonequilibrium 
is  frozen.  Nakamura  (10)  found  a  shift  in  the  distribution  of  some 
radicals  in  a  hydrocarbon  diffusion  flame  with  an  applied  electric  field; 
Place  and  Weinberg  (24)  found  in  a  hydrocarbon  diffusion  flame  that  the 
radiation  intensity  of  OH  increased  with  the  application  of  a  field  while 
the  radiation  intensity  of  CH  decreased. 

In  the  present  investigation  it  cannot  be  determined  if  chemical 
kinetics  are  being  altered  by  the  field  because  the  strong  action  of 
induced  gas  motion  tends  to  overshadow  changes  in  kinetics.  Diffusion 
flames,  however,  are  generally  considered  to  be  reaction  independent 
and  diffusion  controlled,  and,  therefore,  it  is  suspected  that  mechanical 
effects  on  the  flame  are  more  important  than  changes  in  kinetics. 
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The  overall  effect  of  the  field  on  the  combustion  process  can  be 
determined  by  an  analysis  of  the  combustion  products.  Figure  4l  shows 
the  flow  diagram  by  which  a  sample  of  exhaust  gas  was  analyzed  for  CO, 

CQ2,  and  unburnt  hydrocarbon  content.  The  sample  was  withdrawn  from 
the  2-inch  x  2-inch  cross  section  at  the  top  of  the  combustion  chamber 
(see  Fig.  6).  The  tip  of  the  sample  probe,  which  was  a  quartz  tube, 
0.057-inch  O.D.  and  0.025-inch  I.D.,  was  positioned  in  the  center  of 
the  cross  section.  The  gases  were  sufficiently  mixed  at  this  level  that 
the  position  of  the  probe  was  not  important. 

The  CO  and  CO2  concentrations  were  measured  by  infrared  absorption 
analyzers  calibrated  with  a  sample  of  known  concentration  with  nitrogen 
as  a  background  gas.  The  sample  withdrawal  rate  was  less  than  5$  of  the 
total  flow.  Figure  42  shows  that  the  CO2  concentration  for  all  the  flames 
is  practically  independent  of  the  electrode  voltage.  The  CO  concentra¬ 
tion,  Fig.  43,  however,  is  seen  to  decrease  above  voltages  of  around  800 
to  1600  volts,  depending  on  the  flame,  and  then  increase  again  at  volt¬ 
ages  above  2400  volts.  At  the  voltages  of  minimum  CO  concentration,  the 
CO  experiences  a  25$  to  40$  reduction  from  its  no-voltage  concentration. 
The  level  of  concentration  is  so  low  compared  to  the  C02  that  it  cannot 
be  established  whether  the  changes  in  CO  are  reflected  in  changes  in  the 
amount  of  CO2. 

The  unburnt  hydrocarbon  content  of  the  exhaust  was  measured  by  the 
flame  ionization  method.  The  measuring  instrument  was  supplied  with  a 
hydrocarbon  free  mixture  of  hydrogen  and  nitrogen  and  air  which  reacted 
in  a  burner  assembly.  The  sample  was  injected  into  the  burner,  and  if 
hydrocarbons  were  present  in  the  sample,  an  electric  current  was  measured 
in  an  electrode  system  near  the  burner.  The  instrument  could  not  deter¬ 
mine  which  hydrocarbon  was  present,  since  it  is  sensitive  to  both  concen¬ 
tration  and  number  of  carbon  atoms  per  molecule.  For  example,  a  given 
reading  may  correspond  to  a  certain  volume  percent  propane,  or  to  three 
times  as  much  methane.  In  the  present  case,  the  unburnt  hydrocarbon  was 
assumed  to  be  propane,  although  it  is  realized  that  lighter  hydrocarbons 
could  be  present  from  the  products  of  fuel  pyrolysis.  To  remain  within 
the  range  of  the  instrument,  the  sample  rate  was  set  at  around  0.02 $  of 
the  total  exhaust  flow. 

Since  all  the  flames  were  richer  than  stoichiometric,  it  is  not 
surprising  that  unburnt  hydrocarbon  was  present  in  the  exhaust.  Figure 
44  shows  the  measured  concentrations  for  the  flames  investigated.  The 
trends  with  increasing  electrode  voltage  are  similar  to  the  trends  with 
CO  concentration.  A  decrease  is  noted  above  intermediate  values  of 
voltage;  at  high  voltage,  flames  1  and  4  show  a  rise  in  unburnt  hydro¬ 
carbon  following  the  decrease  at  intermediate  voltages. 

The  measured  amounts  of  input  propane  and  exhaust  products  enable 
a  carbon  balance  to  be  made  on  the  combustion  chamber.  Carbon,  which  is 
not  accounted  for  in  the  CO,  CO2,  and  unburnt  hydrocarbon,  is  assumed  to 
be  solid  carbon  formation  and  is  shown  plotted  in  Fig.  45  as  a  function 
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of  electrode  voltage.  The  carbon  balance  indicates  that  solid  carbon  is 
formed  at  zero  voltage,  and  that  the  amount  increases  with  increasing 
voltage.  The  fact  that  there  is  gas  borne  solid  carbon  is  also  supported 
by  experimental  observations;  both  the  gas  sample  filter  and  the  screen 
at  the  exit  plane  of  the  chamber  were  found  to  collect  carbon.  The  in¬ 
terior  walls  of  the  chamber,  although  becoming  black  after  a  time  of 
operation,  did  not  cintinually  build  up  a  deposit  of  carbon.  It  appears 
that  the  following  events  take  place:  propane,  which  normally  at  zero 
voltage  would  escape  unburnt  from  the  chamber,  is  caused  to  pyrolyse  and 
form  free  carbon  as  a  result  of  electrically  induced  mixing  within  the 
chamber.  This  free  carbon  then  mostly  escapes  from  the  chamber,  because 
neither  the  CO  nor  the  COa  concentration  shows  an  increase  with  voltage. 


CHAPTER  10 


EfFECT  OF  THE  FIELD  ON  THE  HEAT  TRANSFER  RATE 


10.1  EXPERIMENTAL  PROCEDURE 

Although  the  great  majority  of  flames  used  in  industry  have  as 
their  purpose  the  supply  of  heat  to  be  beneficially  used  in  some  process, 
the  general  subject  of  flame  heat  transfer  historically  has  not  received 
a  great  deal  of  emphasis  from  either  flame  researchers  or  engineers. 

The  organization  of  the  International  Flame  Research  Foundation  (44)  and 
such  works  as  Hiring* s  (45)  were  therefore  welcomed  by  workers  in  the 
subject.  More  recently  other  works  (e.g.,  46,  47,  48)  concerned  with 
heat  transfer  between  flames  and  enclosing  or  adjacent  surfaces  have 
appeared.  In  the  area  of  special  effects  on  flame  heat  transfer,  sound 
waves  have  been  investigated  and  found  to  be  capable  of  significantly 
increasing  heat  transfer  rates.  Also,  Payne  and  Weinberg  (2)  with  pre¬ 
mixed  flames  in  a  round  tube  showed  that  a  radial  electric  field  could 
increase  the  heat  transfer  rate  up  to  a  factor  of  two. 

In  the  present  investigation,  heat  transfer  rates  were  determined 
by  measuring  the  temperature  rise  and  flow  rate  of  coolant  to  the  elec¬ 
trodes.  The  measurement  system  is  shown  in  Fig.  46  and  the  photographs 
in  Figs.  8  and  9.  Each  electrode  segment,  r.  one-inch-wide  by  one-half¬ 
inch  thick  by  eight-inch-long  brass  bar,  was  drilled  through  longitudi¬ 
nally  to  form  a  l/4-inch  diameter  oil  passage;  a  twisted  ribbon  ol  alumi¬ 
num  was  inserted  into  the  hole  to  promote  mixing  and  the  ends  were  plugged. 
Access  holes  for  the  inlet  and  outlet  flow  of  oil  were  drilled  through 
the  side  of  the  segment  into  which  were  cemented  short  nylon  tubes  which 
served  as  thermocouple  sections  and  tubing  fittings.  The  oil  to  each 
segment  was  supplied  from  a  manifold  containing  a  valve  for  each  segment 
by  which  the  flow  rate  could  be  regulated.  After  passing  through  the 
electrodes,  the  oil  streams  flowed  in  plastic  tubing  to  a  trough  down 
which  the  oil  flowed  into  a  sump.  When  it  was  desired  to  measure  the 
flow  rate  in  any  particular  segment,  the  oil  tube  from  that  segment  was 
manually  diverted  from  the  main  oil  trough  into  a  trough  that  emptied 
into  a  2000  ml  graduated  cylinder.  The  flow  rate  was  then  measured  by 
a  stopwatch. 

The  electrode  segments  when  stacked  up  to  form  the  anode  and  cathode 
were  separated  by  a  l/8-inch-thick  strip  of  asbestos  which  isolated  each 
segment  thermally.  The  stack  of  segments  constituting  either  electrode 
were  held  in  place  by  two  threaded  tension  rods  l/8-inch  in  diameter  en¬ 
cased  in  a  3/l6-inch  O.D.  nylon  tube,  passing  through  a  hole  nr  ur  each 
end  of  each  segment.  The  segments  were  pulled  together  by  tightening 
nuts  on  the  ends  of  the  tension  rods.  The  anode  and  cathode  assembly 


were  then  squeezed  together  against  the  edge  walls  of  -he  combustion 
chamber  by  four  plastic  clamps  on  each  side. 


Two  junctions  of  copper-constantan  thermocouples  were  placed  directly 
in  the  inlet  oil  stream  and  two  junctions  in  the  outlet  oil  stream;  the 
thermal  EMF  which  was  measured  represented  the  temperature  rise  of  the 
oil  stream  as  it  passed  through  the  electrode.  Each  segment  also  had  a 
thermocouple  bead  imbedded  l/4-inch  deep  into  itself  and  cemented  in 
place.  All  thermocouple  leads  went  to  two  switch  boxes  which  were  con¬ 
nected  to  the  measuring  potentiometer  capable  of  reading  to  within  0.0025 
millivolts,  or  about  0.06  degrees  F  of  the  oil  temperature  rise.  Thermo¬ 
couple  wires  from  the  anode  were  at  high  voltage  and  so  were  encased  in 
polyethylene  tubing;  the  anode  thermocouple  switch  box  and  the  potentiom¬ 
eter  were  electrically  insulated  and  provided  with  insulating  extended 
stems  for  switching. 

The  experimental  procedure  for  measuring  the  heat  transfer  rates 
was  as  follows:  When  the  desired  flame  was  obtained  in  the  chamber,  with 
no  voltage,  the  temperature  of  each  electrode  segment  was  read.  It  was 
desired  to  have  the  temperature  of  each  segment  differ  from  its  neighbor 
by  less  than  about  10°F  to  keep  the  heat  transfer  between  segments  belo* 
a  small  value.  The  temperatures  of  the  segments  were  regulated  by  ad¬ 
justing  the  oil  flow  to  the  individual  segments.  It  was  also  desired 
to  maintain  the  temperature  level  of  the  electrodes  at  less  than  about 
50°F  above  ambient  to  keep  the  external  heat  loss  from  becoming  appreci¬ 
able.  The  oil  pump  speed,  bypass  valve,  or  individual  oil  flow  rates 
were  adjusted  to  satisfy  this  temperature  condition.  When  the  observed 
temperatures  no  longer  varied  with  time,  steady  state  was  assumed,  and 
temperatures  were  recorded  for  each  electrode  segment.  A  voltage  was 
then  applied  to  the  anode,  and  the  temperatures  were  measured  again,  etc. 

The  heat  transfer  to  each  segment  was  determined  by  the  equation 


Q  =  p(0.F.)ATCp 


(10.1) 


where  AT  is  the  measured  oil  temperature  rise  in  the  segment,  (O.F.)  is 
the  measured  volumetric  oil  flow  rate,  and  p  and  Cp  are  the  density  and 
specific  heat  of  the  General  Electric  IOC  transformer  oil  coolant  ob¬ 
tained  from  Table  9  below. 

Table  9.  Properties  of  Oil  Coolant 


Temperature , 

°F 

Specific 

Gravity 

Specific  Heat 
BTU  lbm"1  F_1 

70 

0.884 

0.443 

100 

-- 

0.459 

120 

0.866 

__ 

140 

— 

0.481 

The  local  heat  transfer  rate  per  unit  area  was  found  from 


q  =>  -3—  (10.2) 

area 

where  the  area  is  calculated  from  a  segment  height  of  one  inch  and  a 
flame  width  of  5-7/8  inches  and  represents  the  area  of  a  segment  exposed 
to  the  flame.  The  flames  were  found  not  to  extend  completely  across  the 
ideal  width  at  all  levels,  but  since  the  main  study  is  relative  effects 
of  fields,  the  use  of  a  constant  width  is  not  significant. 


10.2  DISCUSSION  OF  DATA 

The  first  task  was  a  determination  of  the  heat  transfer  distribu¬ 
tion  along  the  anode  and  cathode  for  no  applied  voltage,  the  results  of 
which  are  shown  in  Figs.  47  and  48  for  the  four  flames.  Several  state¬ 
ments  can  be  made  regarding  the  results;  first,  the  general  trend  of  the 
heat  transfer  distribution  is  similar  for  all  the  flames.  Second,  the 
heat  transfer  to  the  anode  experiences  a  maximum  at  or  near  the  base  of 
the  flame  and  decreases  upward.  This  is  a  consequence  of  the  flame  being 
nearer  the  anode  at  the  base  and  moving  to  a  more  central  position  in  the 
chamber  at  higher  levels  of  the  flame.  The  trends  also  indicate  that 
radiation  is  secondary  to  convection.  Flame  1  which  is  nonluminous  and 
therefore  would  have  a  relatively  low  emissivity  has  nearly  the  same  heat 
transfer  rates  and  distribution  as  flames  2  and  3  which  are  luminous. 
Also,  the  anode  is  seen  to  have  a  much  higher  heat  transfer  rate  at  lower 
parts  of  the  flame  than  the  cathode.  This  would  not  be  the  case  if  the 
idealized  flame  sheet  had  a  significant  radiation,  since  only  a  small 
fraction  of  the  radiation  would  be  absorbed  by  the  gas  layers. 

The  electric  field  effect  on  the  heat  transfer  is  presented  as  a 
deviation  from  the  heat  transfer  at  no  voltage  and  is  a  more  accurate 
measurement  since  systematic  experimental  errors  are  cancelled  out. 
Figures  49-60  show  the  ratio  of  q  with  voltage  applied  to  q  at  no  volt¬ 
age,  as  a  function  of  electrode  voltage  and  for  several  positions  above 
the  base  of  the  flame.  The  measured  points  correspond  to  the  center  of 
the  electrode  segments.  Flame  3>  "the  richest  flame,  yielded  heat  trans¬ 
fer  measurements  which  frequently  did  not  show  the  same  trends  as  the 
other  flames.  It  is  thought  that  because  of  the  greater  carbon  formation 
in  this  flame  the  effects  of  both  electron  attachment  and  deposition  of 
carbon  on  the  electrode  surfaces  could  be  influencing  the  heat  transfer. 

Since  the  electrically  induced  gas  motion  is  such  as  to  relocate 
the  flame  closer  to  the  cathode  and  farther  from  the  anode,  the  expected 
effect  is  to  decrease  the  anode  heat  transfer  and  increase  the  cathode 
heat  transfer.  Figures  49  and  55  show  that  near  the  base  of  the  flame 
the  anode  heat  transfer  rate  can  be  decreased  up  to  an  85$  reduction  of 
the  heat  transfer  at  no  voltage;  and  the  cathode  heat  transfer  can  be 
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Fig.  47  Heat  Transfer  Distribution  along  Electrodes  for  No  Field,  Flaaes  1  and  2 


Transfer  Distribution  along  Electrodes  for  No  Field,  Flames  3  and  A 
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Plf.  52  The  affact  of  tha  alootrie  field  on  the  heat  transfer 
rate,  j  *  4,125  inches 
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55  The  effect  of  the  electric  field  on  the  heat  transfer 
rate,  y  =  0.75  inches 
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Fig.  59  The  effeot  of  the  eleotrie  field  on  the  heat  transfer 
rate,  y  »  5.25  inohes 
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increased  up  to  a  factor  of  six  within  the  range  of  voltages  applied. 

At  progressively  higher  levels  of  the  flames  the  electric  effect  on  the 
heat  transfer  is  reduced,  as  can  be  seen  from  the  curves  in  Figs.  49-60. 

A  comparison  of  the  results  on  flames  2  and  4,  the  same  mixtures 
but  different  velocities,  shows  a  trend  which  was  observed  before  in 
Chapter  8  on  electric  measurements.  For  several  curves  the  electric 
effect  on  the  higher  velocity  flame  4  is  seen  to  "lag”  that  of  flame  2; 
i.e.,  a  more  rapid  increase  or  decrease  of  heat  transfer  ratio  occurs 
at  a  higher  voltage  for  flame  4.  This  same  effect  was  noticed  in  the 
current  density  curves,  and  it  was  also  noted  from  oscilloscope  traces 
that  flame  4  began  to  flicker  at  a  higher  voltage  than  flame  2.  Thus 
there  is  some  evidence  to  support  the  statement  that  aside  from  induced 
flame  movement,  the  onset  of  electrically  induced  flame  flicker  produces 
a  first  order  effect  on  the  heat  transfer  from  the  flame. 

The  data  on  Figs.  49-60  may  be  plotted  to  show  the  distribution  of 
heat  transfer  ratio  along  the  electrodes  at  constant  voltage,  as  is 
shown  in  Figs.  61-64  for  1200  and  2000  volts.  As  mentioned  earlier, 
both  the  decrease  in  anode  heat  transfer  and  the  increase  in  cathode 
heat  transfer  diminishes  at  higher  levels  of  the  flame.  The  luminous 
flames  also  show  some  values  of  anode  heat  transfer  ratio  greater  than 
one  at  upper  parts  of  the  flame.  This  trend  indicates  that  electron 
attached  carbon  particles  are  active  in  reversing  the  flame  movement 
existing  at  lower  levels  of  the  chamber. 

While  the  maximum  heat  transfer  effects  are  observed  to  occur  at 
lower  parts  of  the  flame,  and  current  densities  are  similarly  a  maximum 
there,  the  current  level  itself  is  not  a  direct  source  of  increased 
transfer.  Electrical  energy  input  to  the  chamber  is  less  than  one  per¬ 
cent  of  the  heat  transfer  measured,  and  the  energy  associated  with  the 
the  work  functions  and  ionization  potentials  is  insignificant. 

At  the  highest  voltage  applied,  the  cathode  heat  transfer  rate  at 
the  first  segment,  Figs.  49,  55  >  is  increasing  with  voltage.  At  higher 
voltages  it  is  therefore  expected  that  even  further  increases  in  heat 
transfer  would  result.  At  these  voltages  the  flame  directly  impinges 
on  the  cathode,  and  of  course  the  model  of  a  flat  flame  is  no  longer 
applicable.  The  flow  situation  more  resembles  stagnation  flow,  as  evi¬ 
denced  by  the  flame  photographs. 


10.3  HEAT  BALANCE 

A  heat  balance  was  calculated  for  each  flame  at  various  electrode 
voltages  and  is  shown  in  Figs.  65  and  66.  The  heat  transfer  to  the 
anode  and  cathode  represents  the  sum  total  of  the  indicated  heat  inputs 
from  Eq.  (10. l).  The  maximum  decrease  in  total  anode  heat  transfer  was 
about  30  to  40$  of  the  anode  heat  transfer  at  zero  voltage.  For  flame  3 
the  maximum  decrease  occurred  at  2000  volts,  and  at  higher  voltages  the 
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anode  heat  transfer  increased,  presumably  because  of  the  effect  of  carbon 
particles  discussed  earlier. 


The  maximum  increase  in  the  cathode  heat  transfer  was  greatest  for 
flames  1  and  2,  about  30  to  40$,  but  a  much  smaller  change  was  found  for 
flames  3  and  4.  In  every  case  the  sum  of  the  total  anode  and  cathode 
heat  transfer  was  reduced  by  the  applied  voltage. 

Because  each  electrode  segment  was  at  a  temperature  higher  than 
ambient,  unavoidable  heat  loss  occurred  through  convection  and  radiation. 
For  each  flame  one  average  value  of  heat  loss  was  determined.  An  average 
electrode  temperature  was  first  determined  by  taking  the  average  segment 
temperature  at  zero  volts,  and  at  the  maximum  voltage  applied.  Radiation 
loss  was  figured  by  assuming  perfect  radiation  from  the  average  tempera¬ 
ture  to  a  sink  at  ambient  temperature.  The  combustion  chamber  was  located 
in  a  fume  hood  with  an  exhaust  fan  drawing  air  over  the  external  surfaces 
at  a  measured  velocity  of  one  foot  per  second.  The  convection  loss  was 
estimated  as  the  heat  loss  from  a  flat  plate  at  the  average  temperature 
in  a  uniform  stream  of  air  at  one  foot  per  second  at  ambient  temperature. 
The  total  external  loss  was  small  but  not  a  negligible  fraction  of  the 
total  heat  transfer. 

The  energy  loss  in  the  exhaust  gases  was  calculated  as  follows: 

First  the  exhaust  temperature  was  obtained  from  twelve  iron-constantan 
thermocouples  spaced  over  the  flow  cross  section  above  the  electrodes 
and  connected  in  parallel  to  a  millivolt  recorder.  The  average  exhaust 
temperature  had  a  rapid  fluctuation  as  indicated  in  the  plot  of  -he  ex¬ 
haust  temperature  versus  voltage,  Figs.  67  and  68.  The  composition  of 
the  exhaust  gas  was  taken  to  be  the  measured  quantities  of  CO,  CO2,  un¬ 
burnt  hydrocarbon,  and  the  input  quantity  of  nitrogen.  All  the  hydrogen 
in  the  burnt  hydrocarbon  was  assumed  to  be  oxidized  to  water  and  exhausted 
as  vapor.  Even  though  some  liquid  water  was  observed  in  the  base  of  the 
chamber,  the  quantity  was  only  a  small  portion  of  the  total  water  formed. 
Solid  carbon  formed  was  assumed  to  be  gas  borne  with  the  exhaust  products. 
The  stack  loss  represents  the  energy  required  to  raise  the  exhaust  prod¬ 
ucts  from  the  inlet  tenperature  to  the  exhaust  temperature,  flowing  at 
the  measured  or  calculated  flow  rate.  The  stack  loss  was  found  to  be 
relatively  constant  with  voltage,  and  an  appreciable  part  of  the  heat 
balance. 

The  chemical  energy  released  by  combustion  was  calculated  from  the 
lower  heating  value  of  the  burnt  propane  and  is  shown  as  a  dashed  line 
in  Figs.  65  and  66.  These  dashed  lines  should  ideally  coincide  with  the 
upper  solid  lines  on  these  figures  in  order  to  account  for  all  the  energy 
transfer.  These  lines  for  each  flame  cross  once  and  show  a  reasonable 
heat  balance  considering  all  the  measured  and  calculated  quantities. 

Flame  3  shows  the  worst  agreement  in  the  total  heat  balance.  The  major 
source  of  error  in  all  the  heat  balances  is  probalby  in  the  input  amount 
of  propane,  which  could  slowly  vary  during  a  test  by  around  one  percent 
of  the  desired  value. 
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10.4  DISCUSSION  OF  RADIATIVE  TRANSFER 

Radiative  heat  transfer  from  flames  depends  to  a  great  extent  on 
the  emissivity.  For  nonluminous  flames  the  energy  is  radiated  in  speci¬ 
fic  bands  of  wavelengths ,  whereas  for  luminous  flames  where  carbon 
particles  are  present  and  achieve  a  high  temperature ,  continuous  radia¬ 
tion  at  a  higher  emissivity  generally  takes  place.  In  some  industrial 
applications  radiative  transfer  is  desirable  and  foreign  particles  are 
injected  into  the  flame  to  increase  the  emissivity  (51).  Although  a 
detailed  study  of  radiative  transfer  was  not  made  in  the  present  study, 
and  there  were  indications  that  radiation  was  secondary  to  convection 
in  this  case,  there  is  in  the  general  case  of  luminous  flames  a  possible 
strong  mechanism  for  an  electric  effect  on  radiative  transfer.  There  is 
evidence  that  a  field  can  greatly  affect  the  growth  and  distribution  of 
carbon  particles  in  flames  (2,  24).  A  theoretical  study  (52)  has  shown 
that  large  variations  in  emissivity  are  obtained  by  changes  in  carbon 
particle  size  and  distribution.  An  indication  of  the  magnitudes  of  the 
effect  of  carbon  particles  on  emissivity  is  given  by  Tables  10  and  11, 
which  show  results  from  Reference  52.  It  appears  that  if  carbon  particle 
size  and  quantity  can  be  effectively  changed  electrically,  then  a  large 
heat  transfer  effect  may  be  expected  in  flames  that  are  primarily  radia¬ 
tive  because  of  carbon. 

Table  10.  Flame  Baissivity  for  Various  Carbon  Particle 
Sizes  and  Wavelength  at  a  Surface  Density  of 
1013  Particles  cm-2 


Particle 

Size 

Wavelength,  microns 

2 

6 

10 

16 

50  A 

0.15 

0.02 

.. 

__ 

200  A 

1 

0.8 

0.4 

0.2 

500  A 

1 

1 

1 

0.96 

Table  11. 

Flame  Emissivity  for  Various  Carbon  Particle 
Sizes  and  Surface  Density  at  a  Wavelength  of 

6  Microns 

Particle 

Flame  Surface  Density, 

Particles  cm"2 

Size 

1011  1012 

1013 

1014  1015 

50  A 

200  A 

500  A 

0.01 

0.02  0.13 

0.2  0.92 

0.02 

0.8 

1 

0.2  0.9 

1  1 

1  1 
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CHAPTER  11 


SUGGESTIONS  FOR  FURTHER  RESEARCH 


The  present  Investigation  was  limited  to  consideration  of  steady 
applied  voltages.  Alternating  electric  fields  could  also  affect  flames , 
and,  because  of  the  field  reversal,  result  in  fundamentally  different 
interactions  with  the  flame.  The  frequency  conceivably  would  be  an  im¬ 
portant  parameter  since  it  would  determine  the  length  of  travel  of  any 
given  ion  before  the  ion  reversed  direction  in  response  to  the  field. 

Three  possible  regimes  of  action  would  occur  when  the  field  and  ion 
mobility  are  such  as  to  reverse  the  direction  of  ion  motion  when  the  ion 
has  traveled  a  distance  less  than  the  flame  thickness,  less  than  the  com¬ 
bustion  chamber  width,  or  greater  than  the  ccmbusion  chamber  width. 

The  high  voltage  equipment  available  could  only  supply  either  steady 
or  60  cycles  per  second  output.  Figure  69  shows  photographs  of  two  flames 
with  60  cycle  voltages  applied  to  the  electrodes.  The  results  were  the 
same  whether  the  fuel  side  electrode  or  the  air  side  electrode  was  im¬ 
pressed  with  the  voltage  and  the  other  side  grounded.  The  sequence  of 
photographs  reveal  a  flame  movement  toward  the  electrode  on  the  air  side 
of  the  flame.  At  intermediate  voltages  one  of  the  flames  possesses  a 
peculiar  wavy  shape. 

Free  carbon  in  the  flame  gases  appears  to  be  significant,  although 
a  flame  movement  toward  the  air  side  was  present  with  an  entirely  non- 
luminous  flame.  The  lower  sequence  of  photographs  in  Fig.  69  shows  a 
flame  relatively  rich  in  fuel  with  a  corresponding  large  amount  of  free 
carbon.  At  the  highest  voltage  the  flame  is  seen  to  be  spread  cut  across 
the  width  of  the  combustion  chamber.  It  is  speculated  that  charged  carbon 
particles  may  have  a  mobility  low  enough  that  they  oscillate  between  the 
electrodes,  inducing  a  turbulence  which  results  in  a  flame  more  character¬ 
istic  of  a  premixed  flame  rather  than  a  diffusion  flame. 

It  is  suggested  that  further  study  be  applied  to  the  interaction 
of  flames  with  controlled  frequency  electric  fields;  in  this  connection 
luminous  hydrocarbon  flames  may  be  particularly  interesting. 
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69  Flame  Photographs  with  Alternating  Voltage 


CHAPTER  12 


SUMMARY  AND  CONCLUSIONS 


Flat  parallel  diffusion  flames  of  propane  and  air  were  formed  in  a 
flat  combustion  chamber;  tbe  wall  on  the  fuel  side  served  as  an  anode 
and  the  wall  on  the  air  side  was  the  cathode  across  which  a  D.C.  voltage 
of  up  to  3200  volts  was  impressed.  Four  flames  were  investigated;  three 
flames  at  equal  velocities  with  various  proportions  of  propane  and  one 
flame  at  a  higher  velocity.  The  leanest  mixturew  was  nonluminous  and 
the  richest  mixture  very  luminous.  The  combustion  chamber  had  segmented 
anode  and  cathode  walls  and  was  instrumented  to  measure  local  current 
density,  local  heat  transfer  rate,  electric  pressure,  exhaust  composi¬ 
tion,  and  exhaust  temperature  as  a  function  of  applied  voltage.  The 
flames  could  be  viewed  through  a  window  in  the  edge  wall  of  the  combus¬ 
tion  chamber. 

The  applied  electric  field  was  found  to  affect  the  flame  through 
three  different  mechanisms.  It  was  shown  that  a  gradient  in  current 
density  along  the  flame  was  responsible  for  an  induced  gas  motion  away 
from  the  flame  in  the  zone  containing  positive  ions.  Such  an  induced 
motion  affected  the  concentrations  of  fuel  and  air  in  the  chamber  so  as 
to  cause  the  flame  to  move  toward  the  cathode. 

Another  interaction  between  the  field  and  the  flame  was  the  onset 
of  a  flickering  condition  which,  with  increasing  voltage,  appeared  at 
upper  parts  of  the  flame  first.  A  parameter  for  describing  flame  insta¬ 
bility  was  set  up  which  Indicated  that  conditions  near  the  saturation 
voltage  of  the  flame  favored  instability.  A  stabilizing  parameter  pro¬ 
portional  to  the  concentration  gradients  in  the  combustion  chamber  w as 
defined  which  indicated  that  the  flame  gets  less  stable  in  the  regions 
further  downstream.  Onset  of  stability  was  detected  by  an  unsteadiness 
in  the  flame  current  observed  on  an  oscilloscope. 

A  third  mechanism  was  the  attachment  of  electrons  to  carbon  particles 
or  molecules  in  upper  parts  of  the  flame.  The  resulting  electrically 
induced  drift  of  the  negative  ions  tended  to  reverse  the  trend  of  induced 
motion  toward  the  cathode  at  lower  parts  of  the  flame. 

The  electric  interactions  with  the  flame  were  found  to  affect  greatly 
the  heat  transfer  rates,  especially  near  the  base  of  the  flame.  At  the 
highest  voltages  applied,  the  flame  was  observed  to  impinge  directly  on 
the  cathode  with  increase  in  local  heat  transfer  over  the  no  voltage  case 
of  up  to  a  factor  of  6.  The  anode  heat  transfer  rate  near  the  base  of 
the  flame  was  found  to  be  reduced  up  to  85%.  At  progressively  higher 
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parts  of  the  flame  the  electric  effect  decreased.  It  is,  therefore, 
concluded  that  heat  transfer  rates  can  he  either  increased  or  decreased 
in  a  controlled  manner  by  the  application  of  fields. 

The  maximum  current  density  from  the  flames  were  only  of  the  order 
of  eight  microamperes  per  square  inch  of  flame  surface  and,  therefore, 
direct  effects  of  current  such  as  recombination  were  extremely  small 
compared  to  the  heat  transfer  rates. 

Appreciable  decreases  in  the  amount  of  carbon  monoxide  and  unburnt 
hydrocarbon  were  found  to  be  obtained  by  the  application  of  voltage; 
solid  carbon  formation  was  found  to  increase.  A  possible  explanation 
for  the  Increase  in  solid  carbon  is  the  electron  attache  *nt  to  carbon 
particles  in  the  pyrolysis  zone  with  a  subsequent  movement  away  from 
the  flame,  favoring  escape  of  the  particle  from  combustion. 

There  were  indications  that  the  heat  transfer  was  convection  con¬ 
trolled,  but  a  possible  strong  mechanism  for  an  electric  effect  on 
luminous  flames  was  outlined.  It  may  be  possible  to  control  the  emissi- 
vity  of  flames  through  the  influence  of  the  electric  field  on  the  forma¬ 
tion,  growth,  or  combustion  of  charged  particles  in  the  flame. 

In  the  present  investigation  an  electric  field  was  applied  across 
the  combustion  chamber  from  the  fuel  side  toward  the  air  side.  From 
the  mechanisms  described,  if  the  applied  voltage  were  reversed  in  polar¬ 
ity,  flame  movement  toward  the  electrode  on  the  fuel  side  would  be  ex¬ 
pected.  In  this  case,  however,  carbon  particles  would  be  present  mostly 
in  a  region  of  positive  ions,  and,  therefore,  ion  attachment  to  particles 
would  add  to  the  effect  of  the  moving  positive  ions. 


APPENDIX 


ONE-DIMENSIONAL  ELECTRIC  EQUATIONS  OF  A  FLAT  FLAME  BETWEEN  ELECTRODES 


Derivation  of  Eg.  (3.1*0 

Integration  of  Eq.  (3.12)  yields 

J+  m  rs  +  constant 

But  J+  «  0  for  s  <  sx,  so  the  constant  above  Is  -rsx,  hence 

J+  ■  r(s-S.L) 

The  current  density  measured  at  the  electrodes  is 

Jo+  "  r(A"si) 

and  also 


JB  -  rA 

Combining  the  above  two  equations 


«i 


which  is  Eq.  (3.14). 


Derivation  of  Eq.  (3.18) 


The  voltage  drop  across  the  ionization  zone  is 


which  is  Eq.  (3.18) 


Derivation  of  Eq.  (3*28) 


From  Eq.  (3.15) 


But 

and 


cLE  n+ 

—  ■  — 

ds  € 
J+  ■  rs 


■ 


it. 

ek+E 


Therefore 


dE  m  rs 
ds  ek+E 


Integrating, 


rs 

V 


+  constant 


But  E  ■  Ef  at  s  •  A,  and  substituting  Jg  ■  rA, 


which  is  Eq.  (3*28). 


Derivation  of  Eq.  (3.29) 

v.  -  j  |V  -  IS  (A2-*2)  ]  d* 

Rearranging, 
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Ef^c 

J*A 

-  1  + 


Evaluating  the  above  equation  results  in 


a*  *  i  .  Jo.*-  fa5“ .  xT 

2  2  l  k€  J  L  JbA  J  L  J_A  J 


which  is  Eq.  (3*29). 


Justification  of  Eq.  (3.35) 

. —  .  ■  —  is  ■  ■■■  ■  ■■  wm  mm  • 

Equation  (3.31*)  can  be  written 


v l^- ]  •  [E'2 - >]* \*T*\ [E'£ * x] - 1]" 


1  p(R-I-A)  ,  E..j 


^  .Sll 


«r 


where  E  is  the  nondlmenaional  field  Strength  E1  ■  E 


For  the  particular  case  when  (R-L)  ■  L  *  50  A  and  E'  ■  2  the  right  side 
of  the  above  equation  is 

50  +  |  +  0.82  +  344  -  | 

The  second  two  terms  are  small  and  therefore  Eq.  (3.35)  is  justified  on 
an  order  of  magnitude  basis. 
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nitrogen,  and  air  were  introduced  separately  at  the  base  of  an  experimental  com¬ 
bustion  chamber  and  burned  in  a  diffusion  flame  sheet  located  between  flat  walls 
of  the  chamber  which  served  as  anode  and  cathode.  The  electrode  walls  were 
instrumented  to  measure  the  local  heat  transfer  rate,  local  current  density, 
and  pressure. 

It  was.fpund  that  the  application  of  a  voltage  difference  across  the 
electrodes  moved  ions  out  of  the  burning  zone  and  resulted  in  a  current  at  the 
electrode  walls.  The  heat  transfer  rate  near  the  base  of  the  flame  was  consider¬ 
ably  increased  on  the  cathode  and  decreased  on  the  anode;  at  positions  further 
from  the  base  of  the  flame  the  electric  effect  was  lessened. 

Flame  distortion  was  thought  to  be  caused  by  electrically  induced  gas 
motion  derived  from  a  gradient  in  current  density  in  the  flow  direction  and 
by  the  onset  of  en  electrically  induced  flame  flickering.  An  analysis  of  the 
products  of  combustion  revealed  that  the  applied  electric  field  acted  to  in¬ 
crease  the  amount  of  unburnt  solid  carbon  and  decrease  the  quantities  of  Tin- 
burnt  fuel  and  carbon  monoxide.  ( 
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